Quantum substitutions of Pisot type, their quantum topological entropy and their use 

for optimal spacing 
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Quantum substitutions of Pisot type and their topological entropy are introduced. 
Their utility as algorithms for optimal spacing is analyzed. 
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II. INTRODUCTION 

The mathematical Theory of Formal Systems, the Theoretical Physics' Theory of Dynamical Systems and the 
Theoretical Computer Science's Theory of Computation may be seen as three different perspectives from which to 
investigate a same object: 

actually, as it has been strongly remarked by John. L. Casti in the 9*'* chapter "The Mystique of Mechanics: 
Computation, Complexity and the Limits to Reason" of one can pass from one perspective to an other through 
the following translation's diagram: 



MATHEMATICS 


THEORETICAL PHYSICS 


THEORETICAL COMPUTER SCIENCE 


formal system 

axioms 
logical inference 
theorems 


dynamical system 
initial conditions 
dynamical evolution 
final state 


computational device 
input 
computation 
output 



that is particularly fruitful as soon as one looks at the phenomenon of Incompleteness. 

Godel's Theorem ^2] stating the existence of propositions undecidable (id est that cannot neither be proved nor 
disproved) in suitable formal systems is seen with indifference, if not with suspicion, by Mathematicians not working 
on foundational issues, the involved undecidable propositions being seen as self- referential intellectual masturbations 
(such as Godel's sentence being a numerical translation, realized through Godel numbering, of a statement asserting 
it own non provability). 

Even more, the same occurs as to Theoretical Physicists whose formalisms so often prefer aesthetical elegance to 
the mathematical rigor, the translation of these formalisms in consistent mathematical terms being seen as a not 
particularly interesting bureaucratic task left to mathematical physicists. 

Contrary, the phenomenon of incompleteness pervades different fields of Mathematics and might have astonishing 
consequences also in Theoretical Physics: 

• as it has been suggested by Geroch and Hartle [1], Markov's Theorem stating the undecidability of the Home- 
omorphism Problem for four-manifolds (see for instance [3|) might compromise our possibility of quantizing 
gravity; 

• as it has been suggested by Agnes and Rasetti Q , the theorems by Boone and Novikov (and their consequences 
such as the theorem stating the recursive undecidability of Dehn's isomorphism problem) stating the unsolvability 
of the Word Problem of suitable finitely presented groups might affect the predictability in suitable classical 
dynamical systems; 

• Berger's Theorem (see 0) stating the undecidability of the Tiling Problem might have consequences in the 
theory of quasicrystals [7[ ; 

• Matyasevich's proof of the undecidability of Hilbert Tenth Problem 0], is a Damocles' sword above the head 
as to all the diophantinc equations occurring in Physics. 

From the other side Thoretical Computer Science has affected Mathematics, allowing to look at the phenomenon 
of Undecidability from a different perspective: 

from the viewpoint of Gregory Chaitin's two Information-Theoretic Incompleteness Theorems [l3| , [HI the existence 
of limitations in the predictive power of formal systems appears quite reasonable and its role in Physics appears natural 
[l^ . in particular as far as Chaos Theory is concerned, as we will show in a specifically dedicated section of this paper. 

It would seem rather natural to think that all this classical stuff could be generalized to the quantum case. 

Unfortunately, the problem of extending Algorithmic Information Theory to the quantum domain is a not yet 
settled issue [H, 0, [13, [11, El- 

In particular the problem, first proposed by Karl Svozil, that almost a decade ago raised our interest in the subject 
during the Phd-studies [11, id est that of formulating and proving quantum analogues of Chaitin's Information- 
theoretic Incompleteness Theorems, is still open. 

Since there exists also a combinatorial approach to Information Theory coarser than the algorithmic approach 
but incomparably easier, it appears natural, as a first step in such a direction, to introduce and study the quantum 
combinatorial information of the quantum analogues of the simplest formal systems, id est those whose inference rules 
are simple substitution rules on the involved symbols [21]. 

In this way we have been naturally led to the first result of this paper, namely the definition of the quantum 
topological entropy for quantum substitutions. 
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Among all the substitutions, furthermore, there exists a particular class, the substitutions of Pisot type, whose 
very useful properties leading to their role in the Theory of quasicrystals and of nonperiodic tilings [13], [11], [13], 
[25I ] are naturally derived from the mathematical properties of a remarkable set of numbers, the Pisot- Vijayaraghavan 
numbers [26j , underlying much of the idolatrous mystique of the golden number t (for this reason we have emphasized 
how the same properties are shared with the poor relative of r, the so called plastic number p.) 

This has led us to the second result of this paper, id est the formulation of a quantum algorithm for optimal spacing. 
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III. TOPOLOGICAL ENTROPY OF SUBSTITUTIONS 



Given a finite alphabet A — \ai, ■ ■ ■ a\A\} (id est a set A such that the cardinaUty \ A\ of A is such that 2 < jAj < Nq) 
let us denote, following [ll|, [IH, and [ISl, by ^+ := UneN+A" the free semi-group generated by A, id est the set 
of all the finite words over A. 

Given x,y G let us denote hy x ■ y the concatenation of x and y, id est the string {xi, ■ ■ ■ ,x\s\,yi, ■ ■ ■ ,y\y\) 
where \x\ denotes the length of the strings x. 

Introduced the set of the sequences over A let us endow A with the discrete topology, and let us endow A~^ 
and with the induced product topology. 

Such a topology over A^+ is the metric topology induced by the following distance: 

d{x,y):^h' , '^1 = ^: (3.1) 

\ ^min{„EN+:x„^„„} , OthCrWlSC. 

In the following we will denote sequences through a bar, id est we will adopt the following notation: 

X := {a;„}„eN+ £ A^'^ (3.2) 
Given a word w — (wi, • • • , Wr) G A'^ and a sequence x € : 
Definition III.l 

w occurs in x at the position n G ; 

Xn — Wi A • • • A Xn+r-l ~ Wr (3-3) 



Definition III.2 

w is a factor of x: 

w <occ X := 3n : w occurs in x at the position n € N+ (3-4) 

In an analogous way, given two words x,y Cz A^ : 
Definition III.3 

X is a factor of y: 

X <occ y ■= 3n e N+ : y„ = Xi A • • • A yn+\x\ = (3.5) 
Given two words x,y €z A'^ let us denote by \x\g the number of occurences of y in x. 

In a similar manner, given a sequence x G and a word y G let us denote by \x\ff the number (eventually 
infinite) of occurences of y in x. 

Definition III.4 

language of length n of x: 

Cn{x) :— {y e A" : y is a factor of a;} (3.6) 

Definition III. 5 

language of x: 

C{x) := \Jnm+Cn{x) (3.7) 

Given a factor y e yl+ of the sequence x G : 
Definition III.6 



frequency of y in x: 



/ (a;) := lim (3.8) 



Let us introduce some remarkable classes of sequences: 
Definition III. 7 

set of the •periodic sequences: 

PERIODIC{A^+) := {x € : (3T e N+ : {xn+T = Xn Vn e N+))} (3.9) 

Definition III.8 

set of the ultimately periodic sequences: 

PERIODICuit{A^+) ■■= {x e : {3T e N+, 3NeN+: (x„+t = a;„ Vn e N+ : n > A^))} (3.10) 
Definition III.9 

set of the Borel-normal sequences over A: 

NORMALiA^^^) := {x G A''+ : /^^(x) = ^ Vy G A+} (3.11) 

Let us now introduce the following: 
Definition III.IO 

combinatorial information function of x: 

Pn{x) = \Cnix)\ (3.12) 

Obviously: 

1 < Pn{x) < ViGAN+,VnGN+ (3.13) 

Proposition III.l 

Properties of the periodic and ultimately periodic sequences: 
1. 

PERIODIC{A^+) c PERIODICuu{A^+) (3.14) 

2. 

3n e N+ : Pn{x) <n ^ xe PERIODICuit{A^+) (3.15) 

Definition III.ll 

topological entropy of x: 

Htop{x) := lim '—^ (3.16) 

where the existence of the limit follows from the subadditivity of the function n log|^| (p„(a;)). 

Remark III.l 
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Clearly the topological entropy of a sequence is different from zero only if Pn{x) grows exponentially with n. 
For instance introduced the following: 

Definition III.12 

Sturmian sequences: 

STURMIAN{A^+) := {x € : p„(x) = n + 1 Vn e N+} (3.17) 
it may be easily proved that: 
Proposition III. 2 

Htopix) = e STURMIAN{A^+) (3.18) 
PROOF: 

Applying the definition IIII . 1 1 1 and the definition IIII. 121 it follows that: 

logi 41 (n + 1) 

Htopix) = lim = Vx e STURMIAN{A^^+) (3.19) 

m 

Remark III. 2 

Let us recall how Andrei Nikolaevic Kolmogorov, in his seminal paper on Algorithmic Information Theory [2^, showed 
with his usual intellectual clearness that there exist three possible approaches to Information Theory: 

• the combinatorial approach in which the information of an object is defined relatively to a context of objects 
distributed uniformly 

• the probabilistic approach in which the information of an object is defined relatively to a context whose objects 
are weighted according to a non necessarily uniform probability distribution 

• the algorithmic approach furnishing an absolute measure of the information of an object, id est not relative to 
any context of different objects of which it is thought to be a member 

As we will now show quantitatively the topological entropy furnishes a measure of the information contained in a 
sequence coarser than the one given by Algorithmic Information Theory. 

From the other side it has the positive feature of being purely of combinatorial nature, id est not to require the 
introduction of concepts from Computability Theory. 

Let us then introduce the following useful partial ordering over A+: 
Definition III. 13 

X is a prefix of y: 

X <py := 3z e : y^x-z (3.20) 

Given a set 5 C A+: 
Definition III.14 

S is prefix- free: 

X -ftp y yx,y e S (3.21) 

Let us now briefly recall some necessary rudiments of Computability Theory [s^] , [3l[ , [13] , [11] . 
Given a partial function C : v4+ A yl+ we can look at it as a total function adding a " non- halting" symbol" (that 
following |ll| we will assume to be the infinity symbol oo) and posing C{x) := oo whether C doesn't halt on x. 
We can then introduce the following: 



Definition III.15 

halting set of C: 

HALTc := {X&A+ ■ C{x) ^ 00} (3.22) 

We will say that: 
Definition III.16 

C is a partial computable function: 

3/ : N A N partial recursive : C{x) ^ string{f [string-^ {x)) (3.23) 

where string : N is the map associating to an integer n the n*'^ word of in lexicographic order and where 

we demand to the mentioned literature as to the definition of a partial recursive function on natural numbers. 
Given a set X C A+: 

Definition III. 17 

X is computable: 

the characteristic function xx of X is computable 
A notion weaker than computability is the following: 
Definition III.18 

X is computably enumerable: 

3C : A+ A A+ partial computable function : X = HALTc (3.24) 

Deflnition III.19 

Chaitin computer: 

a partial recursive function C : ^4+ A such that HALTc is prefix free. 

Definition III.20 

Universal Chaitin computer: 
a Chaitin computer U such that for every Chaitin computer C there exist a constant cu,c € 1R+ such that: 

\/x e HALTc, 3y€A+ : U{y} = C{x) A \y\ < \x\ + cu,c (3.25) 

Given a universal Chaitin computer U and a word x G : 

Definition III.21 

algorithmic information of x with respect to U: 

j^(^) I min{|y| ■.yGA+, U{y) = £}, if {y G : U{y) = x} ^ 0; ^g ^g) 

1 +c», otherwise. 

Given a universal Chaitin computer U and x S A"*": 
Definition III.22 

algorithmic probability of x with respect to U: 

Pu{x) := E (3-27) 

veu-\x) I I 



Definition III.23 
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Chaitin's halting probability with respect to U: 

nu := J2 (3.28) 

In order to show the conceptual importance of Vlu let us recall first of all the following Alan Turing's celebrated 
result [34]: 

Theorem III.l 

Undecidability of the Halting Problem: 

HALTu is not computable 

Chaitin's halting probability codifies the undecidability stated by the theorem IIII.ll in a very useful form owing to 
the following: 

Theorem III. 2 

Chaitin's Theorem on the Halting Problem: 

Given the first n cbits of r^^ i^{Q,u) one can decide whether U{x) < +00 for every x € U^^-^{0, 1}*' 
where ''{0,1} is the nontcrminating symbolic representation with respect to the alphabet {0, 1} of the definition lA. 51 

Let us now introduce the concept of algorithmic randomness (introduced in independent but equivalent ways by 
Martin-L6f, by Solovay and by Chaitin) by defining the algorihmically random sequences as those sequences with 
maximal algorithmic information (and that hence cannot be algorithmically compressed): 

Definition III.24 

algorithmically random sequences on A: 

RANDOM{A^+) := {x G : VC/ universal Chaitin computer 3cc/ E (0, +00) : I{x{n)) > n ~ cu e N+} 

(3.29) 

where x{n) is the prefix of length n of the sequence x. 

Such a notion doesn't depend upon the chosen alphabet in the following sense: 

Theorem III.3 

Independence of algorithmic randomness from the adopted alphabet: 

RAND0M{J^+) = Ta^,A2[RANDOM{J^I+)] VAi, : max{|Ai|, IA2I} G N (3.30) 
where Tai,A2 is the alphabet's transition map of the definition lA. 61 
Definition III.25 

algorithmically random elements of the interval (0,1); 

RANDOM[{Q,l)] := va{RAN DOM{A^+)) , 2 < |A| < (3.31) 
where va is the numerical value map of the definition IA.4I 

Remark III. 3 

Let us remark that, owing to the theorem IIIL31 the definition IIII.25I doesn't depend on the adopted alphabet A. 

The paradigmatic example of a random real is given by Chaitin's halting probabilities: 
Proposition III. 3 
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Halting probabilities are random: 

nu e RANDOM[{0, 1)] VC/ universal Chaitin computer (3.32) 

Let us now observe that clearly: 
Proposition III. 4 

Periodic and ultimately periodic sequences are not random: 

PERIODICM^^^*)<^ RANDOM{A^+) = (3.33) 
PROOF: 

Let a; be a sequence periodic of period T e N+ after a prefix of iV e N digits. 
By the subadditivity of algorithmic information: 

I{x{nT)) < I{x{N)) + I{{xN+i, ■■■ , xt}) + 0(1) Vn G N+ (3.34) 

and is hence bounded. 
■ 

Let us now introduce the following: 
Definition III.26 

recurrent sequences on A: 

RECURRENT{A^+) = {x e : \x\ff = +oo e A+} (3.35) 

Definition III.27 

uniformly recurrent sequences on A: 

RECURRENT^,,{A^^+) :== {x £ RECU RRENT{£'^+) : Vy G 3r e N+ ?7 <occ • • • , a;„+ J Vn e N+} 

(3.36) 

Demanding to [ll| for the proofs we can now state some remarkable property of the algorithmically random 
sequences: 

Theorem III. 4 

Calude-Chitescu Theorem: 

RANDOM{A^+) c RECURRENT{A^+) (3.37) 



Remark III. 4 

The Calude-Chitescu Theorem shows how theorem IIIL 31 has not be equivocated in the following sense. 
Let us consider two finite alphabets Ai and A2 such that Ai C A2; then clearly: 

\x\a = yx e RANDOM{A^I^),ya e A2 - Ai (3.38) 

and hence, according to the theorem 1111.41 

X RANDOM{A^+) yx G RANDOM{A^+) (3.39) 

as it can be appreciated considering the alphabet Ai :— {0,2} and A2 :~ {0,1,2} involved in the construction of 
Cantor's middle third set extensively analyzed in the section [Ul 



Remark III. 5 
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Can the Calude-Chitescu Theorem (id est the theorem 1111.4^ be strengthened by substituting the set 
RECURRENT{A^+) with the set RECURRENTun{A^^) ? 
Up to our knowledge such a question is still open. 

Furthermore, since uniformly recurrent sequences are deeply linked with the notion of almost periodicity, the issue 
is deeply linked with the problem of estimating the algorithmic information's content of an almost periodic function, 
a problem discussed in the section [Cl 

Theorem III. 5 

NORMAL{A^+) c RANDOM{A^+) (3.40) 

Remark III. 6 

The strict inclusion of theorem IIII. 51 implies that though an algorithmically random sequence is always Borel normal, 
there exist Borel normal sequences that are not algorithmically random. 
The classical counterexample is the following: 

Definition III.28 

Champerknowne sequence: 

XChamperknowne '■— 'nGN-^-String^Tl) (3.41) 

where string{n) is the n"* string in lexicographical ordering. 
By construction: 

XChamperknowne & NORMAL{A^^+) (3.42) 

Anyway the Champerknowne sequence is obviously highly algorithmically compressible, being defined through a very 
short algorithm: 

XCharnperknowne ^ RAN DOM (3.43) 

Let us remark that: 
Proposition III. 5 

Link between topological entropy and algorithmic randomness: 
1. 

Htopix) ^ 1 e RANDOM{A^+) (3.44) 

2. 

Htop{x) = 1 ^ X e RANDOM{A^+) (3.45) 
PROOF: 

1. The theorem IIII.4I implies that: 

Pn{x) ^ \A\" Vx e RANDOM{A^*+),yn eN+ (3.46) 

and hence: 

Htopix) = lim - = 1 yx e RANDOM{A^+) (3.47) 

n — >oo 77, 
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2. the fact that the Champernowne sequence of the definition IIIL28I is Borel normal imphes that: 

Pn {xchamperknowne ) = |v4|" VneN+ (3.48) 

and hence: 

fi 

^topi-^Champerknowne) hni — 1 (3.49) 

n^oo 71 

though, as we have already remarked, x is not algorithmically random. 

■ 

Let us now suppose to have a probability distribution f over the alphabet A, id est a map f : Ai-^ [0, 1] such that: 

^/(a) = 1 (3.50) 

and let us introduce the following: 
Definition III.29 

probabilistic information of f: 

Sprobif) :=<-log|^|/>= -^/(a)log|^|/(a) (3.51) 

aeA 

As it is natural for different approaches devoted to formalize from different perspectives the same concept, the 
algorithmic information and the probabilistic information are deeply linked. 

Let us consider, at this purpose, a stochastic process in which, at the n*'' temporal step, a letter x„ from the 
alphabet A is chosen at random according to the distribution f. 

Supposed to have chosen n letters, the random variables xi, • • • ,a;„ are independent and identical distributed so 
that the random vector (xi, • • ■ Xn) has probability distribution: 

n 

/(")(xi,-- - ,x„) := II fix,) (3.52) 
1=1 

Then: 
Theorem III. 6 

Link between the probabilistic information and the algorithmic information: 

hm UBl^ ^ Sif) (3.53) 

where of course: 

</(f(n)) >= J2 f^'^^^Vi^) (3-54) 

Remark III. 7 



A first generalization of the theorem IIII. 61 mav be obtained removing the hypothesis that different letters a^i, • • • , a;„ 
are identically distributed and independent, id est considering the case in which they form an arbitrary stochastic 
process (an arbitrary shift in the language of Abstract Dynamical System Theory [s^, [1^). 

It can be further generalized to arbitrary abstract dynamical systems through a suitable symbolic codification of 
its trajectories. 

The resulting Brudno Theorem [s^], HI], [1^, [13] states that the dynamical entropy of an abstract dynamical 
system , defined as the asymptotic rate of production of probabilistic-information produced by its dynamics, is equal 
to the asymptotic rate of algorithmic information of almost every its orbit. 



Let us now introduce the following basic: 
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Definition III.30 

substitution over A: 

a map a : A\-^ ^+ 

A substitution a over A may be extended to a morphism of ^4+ by concatenation, id est posing: 

a[{xi,--- ,Xn)] := -i^iaixi) (3.55) 
In the same way it may be naturally extended to a map over by posing: 

o-({a;„}„6N) := ■nen<7{xn) (3.56) 

Remark III. 8 

Substitutions are very efficient ways for producing sequences owing to the following: 
Proposition III. 6 

Fixed point of a substitution beginning with a given letter: 
HP: 

aG A 

a substitution such that 17(a) begins with a and \<j{a)\ > 2 

TH: 

3\a{a) e A^+ : a{a{a)) = a{a) 

Given a substitution a it is natural to define its topological entropy as the sum of the topological entropy of its 
fixed point starting with different letters: 

Definition III.31 

topological entropy of a: 

Htop{<T) := Htop{a{a)) (3.57) 

aeA:\a(a)\>2 

Let us now briefly recall some basic notions concerning Matrices' Theory. 
Given n e : n > 2 and a square matrix A £ Mn{^)'- 

Definition III.32 

A is primitive: 

(Ay eNVi,ie {!,••• ,n}) A (3n e N+ : (A")^,- e N+ Vi,i e {1, • • • , n}) (3.58) 

Definition III.33 

A is irreducible: 

$S Hnear subspace of M" : {Av G S W G S) (3.59) 

Let us recall the following basic: 
Theorem III. 7 
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Perron- Frobenius' Theorem: 
HP: 

A irreducible primitive 

TH: 

A admits a strictly positive eigenvalue a (called the leading eigenvalue of A) whose absolute value is greater than 
the absolute value of all the other eigenvalues, a is a simple eigenvalue and there exists an eigenvector with positive 

entries associated to a 

Given a substitution a let us introduce the following: 
Definition III.34 

a is primitive: 

3neN+ : a occurs in cr"(6) Wa,b€ A (3.60) 

Definition III.35 

incidence matrix of a: 

eM|^l(N) := (Af,),, := \a{aj)\a, (3.61) 

Then: 
Proposition III. 7 

(T is primitive <^ is primitive 
Demanding to the appendix [B] for the involved mathematical notions let us introduce the following: 
Definition III.36 

a is of Pisot type: 

Ma- has a leading eigenvalue A such that for every other eigenvalue a one gets A > 1 > |a| 
It may be proved that: 
Proposition III. 8 

Basic properties of substitutions of Pisot type: 
HP: 

(7 substitution of Pisot type 

TH: 

1. CT is primitive 

2. the characteristic polynomial P^ := det{XI — Ma) G Z[A] is irreducible over Q [41] and hence the leading 
eigenvalue a of M^ is a Pisot- Vijayaraghavan number. 

3. in the fixed point a{a) of a starting with a letter a such that |cr(a)| > 2 the frequencies of the letters are given 
by the coordinates of the positive eigenvector associated with the leading eigenvalue normalized in such a way 
that the sum of its coordinates equal 1. 
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where we demand to the appendix |B] for the definition and the remarkable properties of the Pisot-Vijayaraghavan 
numbers. 

AU this stuff may be concretely implemented in Mathematica 5 through the notebook reported in the appendix [El 
Proposition III. 9 

Number theoretic characterization of substitutions of Pisot type: 
HP: 

a primitive substitution 

TH: 

(J is of Pisot type if and only if it its leading eigenvalue is a Pisot-Vijayaraghavan number. 
Example III.l 

The mathematician Leonardo da Pisa, son of the customer inspector Bonaccio and hence later commonly known 
as Fibonacci, considered in his Liber Abaci (published in 1202) the following problem: 
let us suppose to have couples of rabbits such that: 

1. each month a couple of baby rabbits becomes adult 

2. each month a couple of adult rabbits generates a couple of baby rabbits 

Such a situation may be codified representing an adult couple of rabbits with a zero, a baby couple of rabbits with a 
one and introducing the following substitution over the two-letter alphabet {0, 1}: 

a(0) :={0,1} (3.62) 

a{l) :={0} (3.63) 

from which one derives the Fibonacci sequence (j(0): 
{0} 

{0,1} 

{0,1,0} 

{0,1,0,0,1} 

{0,1,0,0,1,0,1,0} 

{0,1,0,0,1,0,1,0,0,1,0,0,1} 

{0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0} 

{0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,0,1} 

{0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,\ 
1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0} 



{0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,\ 

1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,\ 

1,0,0,1,0,1,0,0,1,0,0,1} 
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{0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,\ 
1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,\ 
1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,\ 
1,0,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,0,1,0,1,0,0,1,0,1,0} 



Let us remark that: 

o-"(0) = (t"-1(0) • cr"-2(0) Vn e N : n > 2 (3.64) 

and hence: 

|a"(0)| = Fn+2 VneN (3.65) 

where Fn is the n*'' Fibonacci number (equal to the number of aduh couples of rabbits at the n*'* month) defined 
recursively by 

Fn+2 ■■= Fn+1 + Fn (3.66) 

Fo := (3.67) 

Fi := 1 (3.68) 



Furthermore: 



k"(0)|o = F„+i yneN+ (3.69) 
|cr"(0)|i = F„ Vn e N+ (3.70) 



The incidence matrix of cr is: 



whose leading eigenvalue is the golden number t :— , id est the least Pisot-Vijayaraghavan number of order two 
(see the example IB.1[) . 
Since: 

= ( ^"+1 ^" ) Vn e N (3.72) 

y Fn Fn-1 J 

it follows that a is also primitive and hence is of Pisot type. 

The eigenvector of normalized so that the sum of its components equals to one is {j^, j^)', hence: 

foim) = (3.73) 
1 + r 

= (3.74) 
1 + r 



and hence in particular: 

/o(^(0)) 

Mm) 



= T (3.75) 



that is consistent with the equation 13.691 and the equation 13.701 since: 

Fn+l 



lim = r (3.76) 



Let us remark that since: 

/o(^(0)) 

him) 



T ^ 1 (3.77) 



it follows that the Fibonacci sequence is not Borcl-normal and hence, in particular, is not algorithmically-random, as 
we could have inferred by the fact that it is Sturmian and hence: 



Htop{<j) = Htopim) = (3.78) 



The fact that = r ^ Q implies that the Fibonacci sequence is not ultimately periodic 



Example III. 2 

Let us consider the following substitution over the alphabet of three letters {0, 1, 2}: 

a(0) :={1,2} (3.79) 

a(l) :-{2} (3.80) 

a(2) := {0} (3.81) 
from which one derives the Padovan sequence (t(0) : 
{0} 

{1,2} 

{2,0} 

{0,1,2} 

{1,2,2,0} 

{2,0,0,1,2} 

{0,1,2,1,2,2,0} 

{1,2,2,0,2,0,0,1,2} 

{2,0,0,1,2,0,1,2,1,2,2,0} 

{0,1,2,1,2,2,0,1,2,2,0,2,0,0,1,2} 

{1,2,2,0,2,0,0,1,2,2,0,0,1,2,0,1,2,1,2,2,0} 

{2,0,0,1,2,0,1,2,1,2,2,0,0,1,2,1,2,2,0,1,2,2,0,2,0,0,1,2} 

{0,1,2,1,2,2,0,1,2,2,0,2,0,0,1,2,1,2,2,0,2,0,0,1,2,2,0,0,1,2,0,1,2,1,2,2,0} 
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It may be easily proved that: 

k"(0)| = Pn+2 Vn e N+ (3.82) 
where P„ is the n*'' Padovan number defined by: 

P(0) := P{1) := P{2) := 1 (3.83) 

P„ := P„_2 + P„-3 (3.84) 

Furthermore: 

lim , ''""i")^' ^ hm ^ ^ lim = p (3.85) 

n— >+oo |(T" (0)1 n^+oo P„_|_]^ n^+oo P„_i 

where p := -^^^ — \/69) a ^+ (9+ V69) ^ plastic number, id est the least Pisot-Vijayaraghavan number (see the example 

23 33 

Ea). 

The incidence matrix of a is: 



^ 1 (3.86) 




whose leading eigenvalue is p. 
Since: 

/ Pn-2 Pn-4 Pn-3 \ 

Pn-5 i^«-4 Vn G N : n > 5 (3.87) 

n-l Pn-3 Pn-2 / 

a is also primitive and hence is of Pisot type. 

The eigenvector of M^^ associated to the eigenvalue p and normalized so that the sum of its components equals to 
one is {p ~ 1, "'""t^"^ . and hence: 

' 1+p ' 1+p J 

/o(x(0)) = p-1 (3.88) 

/i(5;(0)) = ^^-^ (3.89) 
1 + p 

/2(5(0)) = (3.90) 
1 + p 

Since: 

/i(x(0)) _ l + p-p2 



/o(x(0)) -l + p2 

/2(S(0)) 1 



/o(x(0)) -l+p2 



^ 1 (3.91) 



^ 1 (3.92) 



it follows that the Padovan sequence is not Borel-normal and hence, in particular, it is not-algorithmically random. 
Since ^ Q and -^j^ll^l^ll = _ij^p-2 ^ Q it follows that the Padovan sequence is not ultimately periodic. 



Example III.3 
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Let us consider the Pell substitution, id est the substitution a over the binary alphabet {0, 1}: 

£7(0) := {0,1} (3.93) 



a(l) := {0,0,1} (3.94) 

from which one derives the Pell sequence: 
{0} 

{0,1} 

{0,1,0,0,1} 

{0,1,0,0,1,0,1,0,1,0,0,1} 

{0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1} 

{0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,\ 
0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1} 

{0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,\ 

0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,\ 

0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,\ 

0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,\ 

0,1,0,1,0,0,1,0,1,0,1,0,0,1} 

{0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,\ 

0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,\ 

0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,\ 

0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,\ 

0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,\ 

1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,\ 

1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,\ 

1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,\ 

0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,\ 

0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,\ 

0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1} 

{0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,\ 

0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,\ 

0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,\ 

0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,\ 

0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,\ 

1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,\ 

1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,\ 

1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,\ 

0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,\ 

0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,\ 

0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,\ 

0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,\ 

1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,\ 

1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,\ 

0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,\ 

0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,\ 

0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,\ 

0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,\ 

1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,\ 

1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,\ 
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1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,\ 
0,0,1,0,1,0,1,0,0,1.0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,\ 

o,i,o,o,i,o,i,o,i,o,o,i,o,i,o,o,i,o,i,o,o,i,o,i,o,i,o,o,i,o,i,o,o,i,o,i,o,i,o,\ 

0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,\ 
0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,1,0,1,0,0,1,0,1,0,0,1,0,\ 
1,0,0,1,0,1,0,1,0,0,1} 



It may be easily proved that: 

|(t"(0)| = an+i VneN+ (3.95) 
where {a„} G are the Pell numbers defined as: 

ao := (3.96) 
ai := 1 (3.97) 
a„ := 2a„_i + a„_2 Vn e N : n > 2 (3.98) 

Furthermore: 

|fT"(0)|o = a„ + a„-i Vn G N+ (3.99) 



The incidence matrix of a is: 



Since: 



|cr"(0)|i = a„ Vn e N+ (3.100) 



= _ (3.101) 




it follows that is a primitive matrix and hence cr is a primitive substitution. 

Since furthermore the leading eigenvalue of Mo- is 1 + \/2 that is a Pisot-Vijayaraghavan number (see the example 
IB.3P it follows that the Pell substitution is of Pisot type. 

The eigenvector of M„ associated to the leading eigenvalue 1 + \/2 and normalized so that the sum of its components 
equals to one is {j^^, 1+^/2 ^' Therefore: 

/o(a(0)) = (3.103) 

/i(a(0)) = (3.104) 
1 + V ^ 

consistently with the equation 13.991 and the equation 13. 1001 since: 

lim = 1 + V2 (3.105) 

n— >+oo a„_i 



Since: 

/o(^(0)) 

him) 

it follows that the Pell sequence is neither Borel normal nor ultimately periodic 



= V2^Q (3.106) 



22 



IV. SUBSTITUTIONS OF PISOT TYPE AS ALGORITHMS FOR OPTIMAL SPACING 

Let us suppose to have to allocate n E N+ : n > 2 petals of a flower so that: 

f . the average distance between neighbors is maximal in order to maximize the exposition of each petal to sun and 
rain. 

2. the distribution of the distances among neighbour petals is the more uniform one. 

Supposed all the petals have unit length so that their vertices belong to the unit circumference of the complex plane 
let us assume that their positions are e**\ • • • e*^" where < 9i < 02 < ■ ■ ■ < Ok < Ok+i < ■ ■ ■ < On < 2tt and let us as 
pose: 

otk ■= dgeodesiciOk+i - Ok) fc=l,---,n (4.1) 

where: 



and where the distance dgeodesic ■ [0, 27r)^ R: 



(4.2) 



ZTT — [Ok+1 — Ok), otherwise. 



is the geodesic distance between e'^'' and 6*^'=+'^ . 

Our optimization problem consists in maximizing the mean of a: 

En 

a := 

n 

while minimizing the variance of a: 



(4.4) 



n 

The solution of such an optimization problem consists in locating the rt*'* petal in the n*'' root of unit: 

e'''^ = rn,k := e^^'= fce{l,-- - ,n} (4.6) 

so that: 

2tt 

Uk = a — — fc = l,---,n (4-7) 

n 

a'^{a) = (4.8) 

To prove it let us observe that, by the eauation l4.3[ any configuration such that max{0fc — 9k-i} < tt and Uk — 2t: 
furnishes the maximum value of the mean of a: 

a = — (4.9) 
n 

Among these configurations the one with minimal variance is, of course, that in which ak = Ok ~ ^fe-i = ^ for every 
k so that the variance of a is equal to zero. 

The roots of the identity have some beautiful properties; first of all: 

Proposition IV. 1 

Cyclotomic identity: 

n 

r„,fc = Vn e N+ : n > 2 (4.10) 

k=l 
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FIG. 1: The S**" roots of unity generate self-similarity. 



PROOF: 



It is sufficient to observe that: 



where we have used the fact that: 



k=l 



k=l 



J2 = ^\ ^ ' Vz e C - {1} , Vn e N+ : n > 2 



fc=i 



1 - z 



(4.11) 



(4.12) 



Let us now denote by polygoun^i the regular polygon having the n*^ roots of unity as vertices. 

Given j, j G {1, • • • , n} such that i j,i ^ j ±1 let us furthermore denote by diagonal{i, j) the segment connecting 

Tn^i Sud TnJ- 

Let us then consider the set of the intersections of distinct diagonals, id est the set of points {diagonal{ii, ji) n 
diagonal{i2,j2) , ii,»2,ji,j2 e {1, • ■ ' ,n},ii ^ ji,i2 ^ j2,H ^ h.ji ^ J2 and let us caU polygo7in,2 the polygon 
having these points as vertices. 

We will say that: 

Definition IV. 1 

the n*"^ roots of unity generate self-similarity: 

polygonn.2 is regular 



Then name of the definition IIV. 1 1 is owed to the fact that if polygonn,2 is regular then there exists a scaling factor 
Tscaiing € (0, 1) and a rotation angle iprotation e [0, 2tt) such that: 



by: 



polygon„,2 = rscaiinge'''^''°^''"''"polygon„^i (4.13) 
Clearly in this case the procedure may be iterated generating a sequence {polygonn.k, k G N+} defined recursively 



polygOnn,k = rscalingt 



^polygonn. 



k-l 



(4.14) 



Example IV. 1 

The 5*'' roots of unity generate self similarity, with r scaling 
and verified looking at the figure [TJ 



l+2r 



and 



:— TT as it can be easily computed 



Let us now suppose to have a countable infinity of petals, so that the uniform distribution given by the equation 

become meaningless. 
Let us then follow another strategy an let us pose the petals in the points: 



PA.fe := exp(i6'A,fc) fc G 



(4.15) 



FIG. 2: The first 10 cusps of Cr- 



where: 

e^.k ■■= Mod2-,[{2TT\f] Vfc e N+ (4.16) 
Let us furthermore call side\^k the segment connecting p^.fe and px,k+i and let us then introduce the following: 
Definition IV.2 

piece-wise linear curve associated to \: 

Cx := UfceN+sirfcfe (4.17) 

Definition IV.3 

oo*'* cusp of C\: 

P\,oo ■= lim p\,k (mod 27r) (4.18) 

Then: 
Tiieorem IV. 1 

Geometrical beauty of Pisot- Vijayaraghavan numbers: 
HP: 



A e PV{A) (4.19) 

TH: 

PA,oo = PA.l (4.20) 
PROOF: 

The proposition IB . 21 implies that if A is chosen to be a Pisot- Vijayaraghavan number, then: 

lim Ox.k = (mod 27r) (4.21) 

k — >+oo 

from which the thesis follows. ■ 

The meaning of the theorem IIV.II is shown in the figure [21 in the figure [3] and in the figure [H 



Remark IV. 1 
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FIG. 3: The first 5 cusps of (7^+^ 




FIG. 4: The first 30 cusps of Cp. 



As for the n*'' roots of unity one can introduce the segment diagonal(i,j) connecting and p\,j for every i,jGN: 
i j,i ^ j ±1 and then considerate the polygon having as vertices the intersections of the diagonals. 
As shown in in the figure O and in the figure [5] no self-similarity is generated in this way. 



As a corollary of the theorem IIV.II it follows that the cusps of Cx are a solution of our optimization problem for a 
countable infinity of petals if and only if A is a Pisot-Vijayaraghavan number. 

Our model is a simplified version of the optimization problem that has led Nature to fix the role of A = t in the 
phyllotaxis of many flowers (see for instance the 11*'* chapter "The golden section and phyllotaxis" of [l^l and (isj). 

Let us now suppose to have a substitution a of Pisot type over the binary alphabet {0, 1} such that |cr(0)| > 2 and 
let us call again A the leading eigenvalue of its incidence matrix. 

Let suppose to choose two fixed angles f3(o):P(i) ^ [0, 27r) such that: 




(4.22) 



FIG. 5: Ct doesn't give rise to self-similarity. 
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FIG. 6: Cp doesn't give rise to self-similarity. 




FIG. 7: The spacing given by the first Ylk^i Fk + 1 — 17711 digits of the Fibonacci sequence, associating spacing r to the letter 
and spacing 1 to the letter 1. 

If we suppose to construct our flower by following the algorithmic procedure: 

^ f Mod^AOk-i + /3(o)), if (^(0))fc = 23) 
\ Mod2.(^fc_i+/3(i)), if (a(0))fe = l ^' ' 

(where ((T(0))fc is the fc*^ digit of the sequence fixed point (t(0) starting with 0) the optimal solution is reached 
asymptotically. 

Examples obtained using the Fibonacci substitution and the Pell substitution are shown, respectively, in the figure 
[7] and in the figure [5] 

Remark IV. 2 

Since to build flowers may be a problem for Nature but is not precisely a concrete problem appearing in common 
life, the optimization problems discussed in this paper might appear as mathematical curiosities with no practical 
application. 




FIG. 8: The spacing given by the first an = 5741 digits of the Pell sequence, associating spacing 1 + y/2 to the letter and 
spacing 1 to the letters 1. 
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The whole matter appears from a different perspective as soon as one reahzes their deep hnk with certain searching 
problems, such as the problem of localizing with the minimum possible uncertainty the minimum of an unknown 
physical observable f whose functional dependence from an other physical observable x is unknown, the only infor- 
mation available being the fact that it is has a unique minimum on an interval (a , b), through a finite number n of 
compounds measurement {xk, fixk)}k=i, each measurement comporting a cost so that the number of measurements 
n has to be minimized too (see the 9*'^ chapter "Optimal spacing and search algorithms" of Ml, the 6*'' chapter 
"Search Techniques and nonhnear programming" of 0|, the I*'' chapter "Basic Concepts" of *46'], the 6*'* chapter 
"Searching" of [47], the 10*'* chapter "Minimization and Maximization of Functions" of 48] and the 10*'' chapter 
"Search and games" of 49]). 
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V. CHAOS AS A PHENOMENON OF UNDECIDABILITY 

Let us start by giving some more information concerning Chaitin's almost mystic fl number. 

With this regard we have, preUminarly, to introduce some basic notion of ComputabiUty Theory over the Reals on 
which, fortunately, agreement has been reached [lOl, [Slj . 
Given a sequence of rational numbers {r„} S Q^: 

Definition V.l 

{r„} is computable: 

3a,6,s : N N recursive : ^ A r„ = Vn G N (5.1) 

b(n) 

Given x e M: 
Definition V.2 

{r„} converges computably to x: 

camp — lim r„ = a; := 3e : N N recursive : {\rk — < 2^" V/c > e(n)) Vrt G N (5.2) 

n — >+oo 

Definition V.3 

X is computable: 

3{r„} e computable : camp— lim r„ = a; (5-3) 

n — ^+oo 

Let us introduce also the weaker notion: 
Definition V.4 

X is computably enumerable: 

3{r„} G increasing , computable : lim r„ = x (5-4) 

n — *H-cxD 

We will denote the set of computable enumerable reals by c.e.(M). 
Then [HI: 

Proposition V.l 

Characterization of Chaitin's halting probabilities: 

{nu U universal Chaitin computer } = RANDOM[{0, 1)] n c.e.(R) (5.5) 

Let us now consider a formal system whose rules of inference form a computably enumerable set of order pairs 
F :~ {(a„, Tn)}neN, the ordered pair (a„, T„) indicating that the theorem r„ is deducible from the axiom a„. 
We will adopt the Mathematical Logic's usual notation: 

Definition V.5 

T is deducible in the formal system F from the axiom a: 

a^pT := (a,T) G F (5.6) 

We can then state the following fundamental (demanding once more to the basic reference [ll| for three different 
proofs): 

Theorem V.l 

Chaitin's V"^ Information-theoretic Incompleteness Theorem: 
HP: 
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F formal system with one axiom a 

ahF"/(a;)>n" ^ I{x) > n (5.7) 



TH: 



3cF e M+ : a hi? "I{x) > n" =^ n < I{a) + cp (5.8) 



Example V.l 

Let us look at a generic substitution a over the finite alphabet A as to a formal system having as axioms the 
letters of A and as inference rules the substitution rule: 

a cr(a) Va G A (5.9) 

X (t{x) yxeA+ (5.10) 

Let us suppose that the substitution a is Sturmian (for instance it might be the Fibonacci substitution of the example 
IIILII) so that the algorithmic information of all the deducible theorems is rather low. 

It follows that may be used only to estimate the information content of objects of rather low information. 

Let us then state also the following: 
Theorem V.2 

Chaitin's 2*"^ Information-theoretic Incompleteness Theorem: 
HP: 



A finite alphabet 

F formal system such that the set Tp of theorems deducible in F is computably enumerable and such that any 
statement of the form "the n*^ cbit of r{o,i}(^;7) is a 0" , "the n*^ cbit of r{o,i}(^;7) is a 1" can be represented in 

Tp and such a statement is a theorem of Tp only if it is true 

TH: 



Tf can enable us to determine the positions and values of at most finitely many scattered cbits of r^Q ]^}(f2[/) 

Let us now analyze how Chaitin Information-theoretic Incompleteness Theorems are at the heart of Chaos Theory. 
Let us, at this purpose, briefiy recall some basic notions. 
Given a finite alphabet A : 2 < |yl| < Nq 

Definition V.6 

Brudno algorithmic entropy of x & -A^+ : 

B{x) := lim Ml^ll (5.11) 

n — >oo fi 

As it has been proved by Brudno himself [s^l : 
Proposition V.2 
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{x e : B{x) > 0} D RANDOM{A^+) (5.12) 

given a classical probability space {X , /i): 

Definition V.7 

endomorphism of {X , /j.): 
T : HALTij, HALTi^, surjective : 

H{A) = fi{T-^A) yA e HALTi, (5.13) 

where HALT/^ is the halting-set of the measure fj,, namely the cr-algebra of subsets of X on which ^ is defined. 

Definition V.8 

classical dynamical system: 
a triple {X , fi, T) such that: 

• (X , /^) is a classical probability space 

• T : HALT^ —t HALT^ is an endomorphism of {X , ji) 
Given a classical dynamical system {X , ji , T): 
Definition V.9 

{X , iJb,T) is ergodic: 

^ n—l 

Zim„^oo-X] '"(^'~'^''('^)) = /"(^)At(^) yA,BGHALTf, (5.14) 
k=o 

Given a classical probability space (X , /x): 
Definition V.IO 

finite measurable partition of {X , fi): 



P = {Po, ■■■ Pn-i}n^n : 

Pi G HALT{ii) i = , • • • , n - 1 

n Pj = <D yi ^ j 
M^- uiL^p,) = 

We will denote the set of all the finite measurable partitions of (X , fi) by V{X , ji). 
Given two partitions P = {Pi}1~o , Q = {Qj}fJo^ e V{X, n): 

Definition V.ll 

P is a coarse- graining of Q {P < Q): 

every atom of P is the union of atoms of Q 

Definition V.12 

coarsest refinement of P = {Pi}^^Q and Q = {Qj}™^^ £ V{ X , n ): 

py Q & r{x,n) 

Py Q := { Pi n i = 0, • • • , n - 1 j = 0, • • • , m - 1} 
Clearly P{X,ii) is closed both under coarsest refinements and under endomorphisms of {X,fi). 
Definition V.13 



(5.15) 



(5.16) 
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Shannon probabilistic entropy of P = {Pi\^Za G 

n-l 

H{P) ■■= -Y.^'(P^)^°SnHP^)] (5.17) 

Given a classical dynamical system CDS := (X , ^ , T) 
Definition V.14 

Kolmogorov-Sinai entropy of CDS: 

hcDS ■■= supp^-p(^x^t.)ltmn^oo^H{VlZoT-''P) (5.18) 

Let us now briefly recall how the orbits of a classical dynamical system may be symbolically codified. 
Let An := {oq, ■ • • , a„_i} be an alphabet of n letters. 
Considered a partition P = {PiYl^Q G 7'(X,^): 

Definition V.15 

symbolic translator of CDS with respect to P: 
: X An-. 

ipp{x) := string{i) : x G Pi (5.19) 
where string{i) is the i"* letter of An in lexicographic order. 
Definition V.16 



k-point symbolic translator of CDS with respect to P: 



:= -^Zli^piT^x) (5.20) 

Definition V.17 



orbit symbolic translator of CDS with respect to P: 
hp 



^p){x) := f^^i^piT^x) (5.21) 

Let us finally introduce the following basic: 



Definition V.18 

CDS is chaotic: 

Then: 
Tiieorem V.3 

Brudno Theorem: 



3PeViX,^) : ijjp'ix) e RAND0M{A"+^) y - h - almost x e X (5.22) 



CDS is chaotic ^ hcos > (5.23) 



Remark V.l 
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Let us remark that, owing essentially to the proposition IV. 21 the fact of having strictly positive Kolmogorov-Sinai 
entropy, though being (according to the theorcm lV.Sp a necessary condition for chaoticity, is not a sufficient condition 
for chaoticity. 

Let us now look at such a dynamical system from a mathematical logic point of view, id est looking at it as a formal 
system. 

Let us suppose that CDS is chaotic and let P E V{X, /i) be a partition such that: 

^'■p°\x) £ RANDOM{A^^^+^) y - fi - almost xe X (5.24) 

Given the initial condition x £ X, the letter tpp^x) G A\p\ is the axiom of such a formal system, 
Let us now observe that the eg nation 15.241 implies that: 

I{il}''p\x)) = On^+oc{n) y — ^ — almost X £ X (5.25) 

Chaitin's 1*'* Information-theoretic Incompleteness Theorem (id est the theorem IV. ip implies that the eventuality 
that the divergence of I{'ip^p\x)) is faster cannot be decided: 

Corollary V.l 

Undecidability underlying chaotic dynamical systems: 

il)p{x) y^CDS I{'ip''p^){x) — On^+cxi{ny^ ^ — fj. — almost X € X (5.26) 



Remark V.2 

Let us remark that such an explanation of the mathematical logic undecidability underlying the phenomenon of chaos, 
holding for discrete-time dynamical systems, is completely independent from Pesin Theorem that, for continuous-time 
d^amical systems, states that the Kolmogorov-Sinai entropy is equal to the sum of the positive Lyapunov exponents 

Also the 2"* Chaitin Information-theoretic Incompleteness Theorem (id est the theorem IV. 2p affects the Theory of 
Dynamical Systems. 

To appreciate this fact it is enough to observe that a universal Chaitin computer U i nay be seen, from the viewpoint 
of the Theory of Abstract Dynamical Systems, as a particular kind of shift (see [H^, [131 and the references therein). 

The Undecidability of the Halting Problem (id est the theorem lIII.2p appears then as a property of such a dynamical 
system ruled by ^lu owing to the theorem IlII. 21 
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VI. QUANTUM TOPOLOGICAL ENTROPY VERSUS QUANTUM ALGORITHMIC RANDOMNESS 

Given the finite alphabet A :— {ai, • • • , a\A\} 
Deflnition VI. 1 

Quantum computational Hilbert space with respect to A: 

Ha ■■= d^l (6.1) 



Given n e N+ : 
Definition VI.2 

71 qu-\A\-its' Hilbert space: 



< a\b> ^ 5a,b Va, 6 e A 



i^«„ := 



\x X x\ 



<x\y>= 5s,g Vx, y e A" 
The quantum analogue of the set A'^ of all strings over A is the following: 
Deflnition VI.3 

Hilbert space of quantum strings with respect to A: 



E 



X >< X 



(6.2) 



(6.3) 



(6.4) 
(6.5) 

(6.6) 



(6.7) 
(6.8) 



< x\y > ^ Ss^g Vx, y e A+ 
The quantum analogue of the set A^+ of all sequences over A is the following: 
Deflnition VI. 4 

Hilbert space of quantum sequences with respect to A: 



n 



Ha 



/ dxlx >< x\ 
Ja''+ 



(6.9) 



(6.10) 



(6.11) 



< x\y > = S{x — y) Vx, y G A^ 



(6.12) 



where dx is the Lebesgue measure over A^+ and where the S in the right hand side of l6.12l is the Dirac delta tempered 
distribution (defined in the appendix [X)) . 

Given a quantum sequence |'0 >€ Ti^^ the natural quantum analogue of the definition IIII.lOl is the following: 
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Definition VI. 5 

combinatorial quantum information function of \ip >: 



p„(|V >) / dx\ < > (6.13) 



Remark VI. 1 



The definition IVI.5I is nothing but the expectation value of the classical combinatorial information function when a 
measurement of the sequence operator a: ^ is performed. 

The natural quantum analogue of the definition IIII.llI is then the following: 

Definition VI.6 

topological quantum entropy of \^ >: 

H,.M » := lim (6.15) 

Clearly: 
Proposition VI. 1 

Htop{\x >) = Htopix) V.T e (6.16) 
PROOF: 

It is sufficient to observe that: 



p„{\x>) = / dy\<y\x>\^pn{\y>) ^ / dyS^ {x - y)pn{\y >) ^ Pn{x) (6.17) 



Example VI. 1 

Given the binary alphabet {0, 1} let us consider the quantum sequence: 



Then: 



Since obviously: 



|V'>:= ^(10000000000- ••> +11111111111 •••>) (6.18) 
v2 



,. , , p„(|0000000000--->)+p„(|llllllllll--->) 
PnW >) = n (6-19) 



/:„(oooooooooo---) = {•■LiO} (6.20) 

p„(0000000000---) ^ 1 (6.21) 
/:„(1111111111---) = {-^^il} (6.22) 



^ where obviously the sequence operator is defined by: 

S|x > := x\x > (6.14) 
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p„(llllllllll---) = 1 (6.23) 

it follows that: 

Pnili^ >) = 1 (6.24) 

and hence: 

Htopil^P >) = (6.25) 

Example VI. 2 

Let us consider a gaussian wave-packet centered in the sequence 1 0000000000 • • • >: 

|<x|V'>P = exp(-d2(^^ 0000000000 •••) (6.26) 

where d is the metric defined by the equation l3.1l 
Then: 

p„(|?/,>) ^ / dxexp{-(f{x,\QOOOOOOOOO--->)pn{\x>) = V / dxexp{2-^''+^'>^)pn{\x >) 

Ja'''+ fr[J{seA"+:x(k+i)={-'^^,o)-i} 

(6.27) 

Let us now introduce the following: 
Definition VI. 7 

Coleman- Lesniewski operator: 

tluANDOM ■= dxXj,^^ooMiA«.)\^ >< ^1 (6-28) 

where: 



Xsix) := I I' '[l""^: (6.29) 
I 0, otherwise. 

is the characteristic function of a set S. 

It may be easily verified that YIrandom ^ projection operator. It appears, then, natural to define the space of 
the algorithmically random quantum sequences as the subspace on which the Coleman-Lesniewski projects: 

Definition VI.8 

subspace of the algorithmically random quantum sequences with respect to A: 

RANDOMinT^) := t] = {|V >€ : TT |V >= |V >} (6-30) 

^ ^ ' i-i-RANDOM ^ ^' ^ i-i-RANDOM' if J \ ^ 

We can then prove the quantum analogue of the proposition IIII.5I 
Proposition VI. 2 

Link between quantum topological entropy and quantum algorithmic randomness: 
1. 

fftopdV' >) 7^ VIV' >e RANDOM{n®/+) (6.31) 

2. 

Htop{\i^ >) ^ ^ I^A >e RANDOM{nT+) (6.32) 
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PROOF: 

1. Clearly: 

Pn{\lp >) = / dx\ < x\tlj > \^Pn{\x >) = 

J RANDOM{A-'+ ) 

[ dx\ < xlip > \^\A\" > \f\ip >e RANDOM {'Hf^+) (6.33) 

J RANDOM{A"+ ) 

2. the quantum Champerknowne sequence is not algorithmically random: 

\xChamperknon,ne > ^ RAN DOM (nT^) (6.34) 

Anyway: 

^topi]'^ Champerknowne ^) — H fQp{x (j hamper know7ie) — 1 (6.35) 



Remark VI. 2 

The proposition IVI. 21 shows that the considerations of the remark BlI . 2 1 may be thoroughly extended to the quantum 
case: 

the quantum topological entropy furnishes a measure of the quantum information contained in a quantum sequence 
coarser than the one given by Quantum Algorithmic Information Theory [l^, [3, [III, [13, [3, [isj - 

Since no general a gree ment has been reached in the scientific community as to the right quantum analogue of the 
definition IIII.21I (see |4C)j for a discussion of the involved issues), the combinatorial approach has the advantages not 
to involve Quantum Computability Theory with all its still not yet settled issues. 

Example VI.3 

Given fii, 0,2 G RAN DOM[{0, l)]nc.e.(R) and a finite alphabet A let us consider the quantum sequence of qu-|A|-its: 

IV' > ^(|rA(r!i) > +|r^(r!2) >) (6.36) 

where ta is the nonterminating symbolic representation with respect to A of the definition IA.5I 
Obviously: 

IV >e RANDOM{7f/) (6.37) 
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VII. QUANTUM ALGORITHMS OF SUBSTITUTION 

Given two finite alphabet A and B and a map / : —>■ there exists a canonical way of constructing a quantum 
algorithm associated to f as the linear operator / : ^ Ti.^^ acting as f on the computational basis (about 

some remarkable considerations concerning the implementation of such an algorithm see the section 4.9.13 "Quantum 
implementations" of [ssj). 

Applied in particular to a substitution a such a strategy leads to the following: 

Definition VII. 1 

quantum algorithm of the first kind associated to a: 
the linear operator over U 7i^^^ defined by the following action on the computational bases: 



a\x > := |(t(x) > f g ^+ 
a\x > := \a{x) > x e 

We will refer to cr as to a quantum algorithm of substitution or, more concisely, as to a quantum substitution. 
Clearly: 

Proposition VII. 1 

1. 

a\x > = (»f2^\a{xi) > VfeA+ 

2. 

a\x > = (E>nefi+W{xn) > Vi e 



(7.1) 
(7.2) 



(7.3) 



(7.4) 



3. 



4. 



a is not self-adjoint 



(7.5) 



xi£A+ 

PROOF: 

1. It is sufficient to observe that: 

a|f>= \aix)>^ |a(-l^ix,)>= \ -f^^ a{x,) > = ®gi|a(x,)> ^x e A+ 

2. In an analogous way: 

a\x > = \a{x) > = \a{-„(zfiXn) > = \ ■nefi+ (^ixn) > = (^nenJcr^Xn) > G A^^ 

3. It is sufficient to observe that: 



(7.6) 



(7.7) 



4. Clearly: 



x£A+ y^A+ X&A+ y&A+ 



from which the thesis follows. 
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Let us now suppose that a G ^ is such that |cr(a)| > 2. Denoted as usual with (T(a) the fixed point of of cr starting 
with a let us observe that: 

Proposition VII. 2 

a\a{a) > = \a{a) > (7.10) 
PROOF: 

It is sufficient to observe that: 

a\a{a) > = \cr[a{a)] > = \a{a) > (7.11) 

■ 

Let us introduce also the following: 
Definition VII.2 

symmetric state of n qu-\A\-its: 

The Mathematica 5 notebook of the section |E] may be used to compute the action of quantum substitutions on 
symmetric states. 

Example VII. 1 

Let a denote the Fibonacci substitution of the example IIII.ll Then: 

a\S, 1 > = a[-^(|0 > +|1 >)] = -i=(|01 > +|0 >) (7.13) 

a^\S, 1 > = ^(|010 > + |01 >) (7.14) 
v2 

<T^|S', 1 >= ^(101001 > +|010 >) (7.15) 
v2 



and in general: 



and hence: 



a"\S, 1 > = ^(|(t"(0) > + \<7''-\0) >) Vn e N+ (7.16) 
v2 



lim ct"|S',1>= ^(|a(0) > +|a(0) >) = V2\a{0) > (7.17) 



Furthermore: 



&\S,2> = (7^(100 > +|01 > +|10 > +|11 >) = ^(10101 > +|010 > +|001 > +|00 >) (7.18) 
a^\S, 2 > = i(|010010 > +|01001 > +|01010 > +|0101 >) (7.19) 
a^\S,2>= i(|0100101001 > +101001010 > +101001001 > +1010010 >) (7.20) 
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Remark VII. 1 

The quantum algorithms of substitution of the first kind have the great disadvantage of loosing the connection with 
the powerful and beautiful Mathematics of Pisot-Vijayaraghavan numbers underlying the incidence matrix and are of 
doubtful practical utility. 

This suggests the introduction of quantum algorithms of substitution of the second kind. 

Given a substitution a : 
Definition VII.3 

quantum algorithm of the second kind associated to a: 

the linear operator a represented by M„ in the computational basis of li-A- 

We will refer again to ct as to a quantum algorithm of substitution or, more concisely, as to a quantum substitution. 
Then: 

Proposition VII.3 

HP: 



G Pisot substitution with leading eigenvalue A 

a&A: \a{a)\ > 2 
a quantum algorithm of the second kind associated to a 



(7.21) 



TH: 



where: 



lim (t" = |eA >< e^l 



a\ex > = X\ex > 

< ex\ex > = 1 
PROOF: 



(7.22) 

(7.23) 

(7.24) 



Since: 



a = X\ex >< ex\ + ^ y\ey >< Cy 

yeConiX) 



(7.25) 



it follows that: 



Since A e PV{A): 



yeCon{\) 



lim y"\ey >< e„| = 



(7.26) 



(7.27) 



from which the thesis follows. 
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Example VII.2 

The quantum algorithm of the second kind associated to the Fibonacci substitution is clearly the operator a : 7^{o,i} ^ 
H{o,i} such that: 

ct|0 > = |0 > +|1 > (7.28) 

a\l > = |0 > (7.29) 
From the example IIII.ll we may immediately infer that: 

CT = r|er >< e^l + (--)|e^ >< e^l (7.30) 



(where we have used the fact that = r + 1) so that the probability that a measurement of the qubit operator in 
the state limn^+oo o'"|0 >= \er > gives as result zero is: 

Prum„^^^a-\o>iO) = (7.33) 
while the probability that such a measurement gives as result one is: 

i'nim_+„*"|o>(l) = ^ (7.34) 

Example VII.3 

The quantum algorithm of second kind associated to the Padovan substitution a is clearly the linear operator 
: "^{0,1,2} '-^ "^{0,1,2} such that: 

d-|0 > := |1 > +|2 > (7.35) 

a\l > := |2 > (7.36) 

a\2 > ■= |0 > (7.37) 
From the example IIII.ll we may immediately infer that: 

a\ep > = p\ep > (7.38) 

' ^ > + , ^ ^ 1 > + , 2 > 



V(p2-l)2 + (l+p-p2)2 + l' ^(p2_l)2 + (l+p_p2)2 + l' ^(p2 _ 1)2 + (1 + p _ ^2)2 + ^ 

(7.39) 

(where we have used the fact that = p + 1) so that the probability distribution of the outcome of a measurement 
of the qutrit operator in the state lim„^+oo ct'^IO > = |ep > is: 

Pni™„_,^."|o>(0) = n2i^n'/r .2,2^1 (7-40) 



(p2- 


1)2 4 


-(1 + 


p-p2)2 + l 




(1- 


f p- 


p2)2 


(p2- 


1)2 4 


-(1 + 


p-p2)2 + l 






1 





^'nim„_+^a"|o>(l) = j-2 — 1^2 — ;^WVT (7-41) 



Pn..„.,..^|0>(2) = (,2_l)2 + (l+,_,2)2 + l (7-42) 
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Example VII.4 

The quantum algorithm of the second kind associated to the PeU substitution is clearly the operator a : Ti^Q ij t—i- 
Tijo.i} such that: 

<t|0 > := |0 > +|1 > (7.43) 

a\l > := 2|0 > +|1 > (7.44) 
From the example IIII.3I we immediately infer that: 

a = (1 + V2)\e^^^ >< e^+^l + (1 - V2)\e,_^ X e,_^\ (7.45) 



|ei_v^>- -\/olO>+ — 11> (7.47) 



so that the probability distribution of the outcome of a measurement of the qubit operator in the state lim„_>+oo <5'"|0 > 
= |ei+V2 > 

Pr-u„,_^^*.|o>(0) = ^ (7.48) 
^'nim_+^*"|o>(l) = I (7.49) 



Remark VII. 2 

Quantum algorithms mimicking the recursive structure of Fibonacci numbers have been introduced in . In our 
modest opinion it doesn't seem, anyway, to be any connection between those remarkable algorithms and the quantum 
substitutions of Pisot type discussed in this paper. 
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Remark VII.3 

The idea of taking into account suitable (specifically almost periodic Schrodinger) operators associated to substitutions 
has been already considered in the Mathematical-Physics' literature concerning quasicrystals (see for instance [60| . 
[HH and references therein). 

In order to allow the reader to appreciate the possible conceptual connections a brief review of the basic notions is 
presented in the section [Dl 

From the mathematical side the great achievement obtained by Barry Simon and coworkers, namely the discovery 
that almost periodic Schrodinger operators has commonly continuous spectral measures supported on Cantor sets |62l | 
has never been looked from an information theoretic perspective. 

As we have already observed in the remark Illl.Si this would be a necessary step in order to generalize the Calude- 
Chitescu Theorem (id est the theorem IIII.4P clarifying the algorithmic information theoretic status of uniformly 
recurrent sequences. 

The attempt to make some step forward in this direction is presented in the section [Cl 

Remembering, furthermore, that an almost periodic function / : M.^ i-^ C can be seen as an ergodic stationary 
random function over the probability space {T, Hull{r), fiHaar(r)), T :— Tgf being the closure of all the shifts of f 
(where of course T^f{y) := f{x + y)) (see the section 7.2 "Periodic and Almost Periodic Potentials" of [il]) such a 
task could be part of a more ambitious ompletely unexplored research project investigating algorithmically-random 
Schrodinger operators that, anyway, is far beyond the purposes of this paper. 

The natural analogue of the definition IIII.31I for the quantum algorithm of first kind associated to a substitution is 
the following: 

Definition VII.4 

quantum topological entropy of a: 

Htopicr) ^^AHtop{\^ >) , IV' >e : < a{a)\^P >^ 0} (7.50) 

aeA : |cr(a)|>2 



Remark VII.4 



Quantum analogues of the theorem IIII.6I expressing the link between the quantum probabilistic information (id est 
Von Neumann entropy) and the quantum algorithmic information have been established in the same papers in which 
this latter notion has been introduced [11], [ll], [11], [13], [ll], [ll]. 

After having defined chaotic a quantum dynamical system whose orbits, simbolically codified in a suitable way, 
are algorithmically-random, as we have observed in the remark [IIL3 the next conceptual step would then consists in 
its generalization to arbitrary quantum stochastic processes and then, through a suitable symbolic codification, to 
arbitrary quantum dynamical systems resulting in a quantum analogue of the Brudno Theorem (id est the theorem 
IV. 3p stating that the vanishing of the quantum dynamical entropy of a quantum dynamical system (defined as the 
asymptotic rate of production of probabilistic information of such a dynamical system) is a necessary condition for 
chaoticity also at the quantum level. 

The first problem, with this regard, consists in the fact that, according to whether or not one assumes that 
measurements are performed on the quantum dynamical system during its dynamical evolution, one results in different 
notions of quantum dynamical entropy: respectively the Connes-Narnhofer-Thirring entropy [g^ ]. [gsI ]. [66l ]. [btI ] and 
the Alicki-Lindblad-Fannes entropy [68| . 

A remarkable result going in the right direction has been obtained in [69] though, in our modest opinion, it is not 
yet a Quantum Brudno Theorem, since it in no way takes into account individual trajectories and their eventual 
quantum algorithmic randomness that, according to the analysis presented in the section |Vl is the key point as to the 
underlying undecidability phenomenon. 

The definition IVII.41 involving only the combinatorial approach to quantum information, might be a tool useful in 
order to proceed in the indicated direction. 
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VIII. A QUANTUM ALGORITHM FOR OPTIMAL SPACING 

In the section IIVI we have seen that the optimization problem consisting in aUocating n e N+ : n > 2 petals of a 
flower so that: 

1 . the average distance between neighbors is maximal in order to maximize the exposition of each petal to sun and 
rain. 

2. the distribution of the distances among neighbour petals is the more uniform one. 
has as solutions the n*''-roots of unity: 

r„,fc := e*^'^ k = l,--- ,n (8.1) 
Given a generic finite alphabet A and n e N+ [70j : 
Definition VIII. 1 

quantum Fourier transform on n qu-\A\-its: 

the linear map F over 7i^" defined by its matrix's elements in the computational basis: 

< y\F\x > := -J=riA\^^k„^^s)k„^^g^ Va;,y e A",Vn G N+ (8.2) 

where koc denotes the local lexicographic number on A" defined as: 

hoc{-l^iai) = (8.3) 



hocrk=ia\Ai) = \Ar-l (8.4) 

Let us face, anyway, the more interesting optimization problem in which the number of petals is countable. 
Let us now suppose to have a substitution cr of Pisot type over the binary alphabet {0, 1} such that |<t(0)| > 2 and 
let us call again A the leading eigenvalue of its incidence matrix. 
Let suppose to choose two fixed angles /3(o), G [0, 2tt) such that: 

Mod^A^] = A (8.5) 

P{i) 

Let us then construct our flower according to the following algorithmic procedure: 

1. set n=0 

2. label[start] 

3. prepare the state |0 > on the 1*^ register 

4. apply n times the operator & on the I*'' register 

5. copy the status of the 1*'* register on the 2*'' register 

6. perform on the 2*'' register a measurement of the qubit operator 

7. pose the rt"* petal at the angle: 

^ _ I Mod2Tr {6,1-1 + /3(o)), if the result of the measurement was 
1 Mod2-K{0n-i + /3(i)), if the result of the measurement was I 
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8. increase n of one unit 

9. goto the label [start] 

It is important to stress that the fifth step of such an algorithm doesn't violate the No Cloning Theorem as one 
can infer taking into account the precise formulation of such a theorem : 

Theorem VIII. 1 

No Cloning Thorem: 
HP: 

H Hilbert space 

TH: 

■.Hi^'H^H unitary : JJ\ipi >= \tpi > > AC/|V'2 >= > «'|V'2 > A >^ |V'2 > A < tpi\^2 >y^ 

(8.7) 

PROOF: 

Let us suppose ad absurdum that such a unitary operator exists. Then: 

< iJl\Lf^U\^2 >={< V'll'^ < V'l|)(l^2 > ®|V'2 >) = (< V'l|^2 >)' = < V'llV'2 > (8.8) 

and hence: 

< ^1^2 > - (8.9) 

contradicting the hypothesis that the states {ipi > and |?/'2 > are non-orthogonal. ■ 

Theorem I VIII. 1 1 forbids the simultaneous cloning of non-orthogonal states, but it doesn't forbids the simultaneous 
cloning of distinct orthogonal states. 

Given a finite alphabet A and considered the cloning function f clone '■ ^ one can immediately introduce 
the quantum algorithm fdone '■ Ti^"*" i— > 7i^^ acting as fdone on the computational basis: 

fclone\x>:=\x-X> X & A+ (8.10) 

perfectly compatible with the theorem lVIII.il 

The presented algorithm reaches asymptotically the optimal solution. 
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APPENDIX A: GENERALIZED FUNCTIONS ON THE SPACE OF ALL SEQUENCES OVER A FINITE 

ALPHABET 

Given a topological space {X,T) fn\ and a subset of its S C X: 
Definition A.l 

set of the limit points of S: 

LP{S) := {xeX : {O e T A x e O) ^ O n S ^ 9} (Al) 

Definition A. 2 

S is a Cantor set: 

{S = S) A {S° = 0) A {LP{S) = S) 

where S is the closure of S while S° is its interior. 

Given a finite alphabet A = Wi, ■ ■ ■ ct\A\} (id est a set A such that the cardinality of A \A\ g N+) we will denote, 
following the notation of [IH, [131 and jll], by := U„gN+^" the free semi-group generated by A, id est the 
set of all the finite words over A. 

Given x,y G ^4+ let us denote hy x ■ y the concatenation of x and y, id est the string {xi, ■ ■ ■ ,a;|j|,yi, • • • ,y\y\) 
where \x\ denotes the length of the strings x. 

Introduced the set j4^+ of the sequences over A let us endow A with the discrete topology, and let us endow A'^ 
and A^^+ with the product topology. 

Such a topology over A^+ is the metric topology induced by the following distance: 

dix,y) = l^' , '^f = ^? (A2) 

[ ^min{„EM^:x„^y„} , OthCrWlSe. 

The Borel-CT algebra ) associated to such a topology is the cr-algebra generated by the cylinder sets: 

Cs ■■= {y e A'*+ : y{\x\) = x} x € A+ (A3) 

where y(n) denotes the prefix of length n of the sequence y. 
We can then introduce the following: 

Definition A. 3 

Lebesgue measure associated to A: 

the probability measure ULebesgucA over the measurable space (A^+,S(A^+)) such that: 

IJ-Lebesgue,A{Cs) ■= |^||-| (A4) 

Given / e L\A^^+,d^iLe besgucA) pose dx '.= diiLebesgue,A and hence: 

dxf{x) := / dflLebesgueA{^)f{x) (A5) 

A"+ 

The applicability of Probability Theory to the real world lies, from a foundational perspective, on the following: 
Theorem A.l 

On the Foundation of Probability Theory: 

fiLebegue,A[R^NDOM{A''+)] = 1 VA : 2 < |.4| < Ho (A6) 

Given a sequence x £ A^+ : 
Definition A. 4 
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numerical value of x: 

va{x) := l^-rjuT- (A7) 

n=l ' ' 

Remark A.l 

Let us remark that the map v : A^'^+ t-^ [0, 1] is surjective but it is not injective since: 
Proposition A.l 

va{x ■ a^ ■'^^\g\^2 a\A\) = aj+i -^1^1+2 ai) Vf € A+, Vi{l, • • • , |A| - 1} (A8) 

Example A.l 

Let us suppose that A :— {0, 1, 2, 3, 4, 5, 6, 7, 8, 9}. Then we are taught from childhood that: 

0.9 = 1 (A9) 



0.19 = 0.2 (AlO) 
0.29 = 0.3 (All) 



and so on. 

Let us then introduce the following: 

Definition A. 5 

nonterminating symbolic representation with respect to A: 
the map : [0, 1] ^ A^+ : 

lex{xn)-. 



^"^^1^^ |^|« ) -nGN+Xn (Al2) 



with the nonterminating choice in the cases of the proposition lA.ll 

Given two finite alphabets Ai and A2 we can then introduce the following: 

Definition A. 6 



alphabet's transition map from Ai to A2: 



the map Tai,A2 ■ ^i*^ ^2^- 



The structure of topological space that we have given to ^4'^+ is sufficient, as it is well known, to define limits and 
derivatives. 

Given a smooth (id est infinitely differentiable) map / : 1— > C and n, m G N: 
Definition A. 7 

|1/|1„„„ := sup \{vA{x)r — f{vA{x))\ (A14) 

Definition A. 8 

Schwartz space of rapid decrease functions over A^+ : 

5(A^+) := {/ -.A^^+i-^C smooth : ||/||„,™ < +cx)Vn,TO G N} (A15) 
The family of seminorms ||/||„.m induces a natural topology over 5(A^+) that can be used to introduce the following: 



47 



Definition A. 9 

space of tempered distributions over : 

S' {A^+ ) := the topological dual of S{A^+ ) . 

Given S e A^+: 
Definition A. 10 

Dirac delta in x: 

5x := the linear functional := f{x) 

It may be easily proved that: 
Proposition A. 2 

5s,eS'{A^+) Vse (A16) 
In this paper we will adopt the usual notation: 

/ dxf{x)5{x-y) := 5y{f) (A17) 
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APPENDIX B: PISOT-VIJAYARAGHAVAN NUMBERS 

Given x € C: 
Deflnition B.l 

X is algebraic: 

3P e V[Q] : Fix) ^ (Bl) 

where we have denoted with V[K] the hnear space of all the monic polynomials with coefficients on the generic ring 
K. 

We will denote the set of all the algebraic numbers by A. 
Definition B.2 

X is algebraic integer: 

3P e V[Z] : P{x) = (B2) 
We will denote the set of all the algebraic numbers by A^. 

Given n G N+ let us denote with Vn [IK] the set of all the monic polynomials of degree n with coefficients belonging 
to the generic ring K. 
Given x E A: 

Definition B.3 

degree of x: 

deg{x) min{n G N+ ; P{x) = P G VnM} (B3) 

Definition B.4 

minimal polynomial of x: 

Px— minimal ^ '^de(;(2;) [^] ■ Px—minimal{-^^ (-^^) 

Definition B.5 

set of the conjugates of x: 

Con{x) :^ {y eC : y ^x A Px-minimai{y) = 0} (B5) 
We can then state the following [z2| , [z3| : 
Proposition B.l 

Fundamental property of algebraic integers: 
HP: 



X G 



TH: 
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+ y" G Z Vn G N+ : n > deg{x) 

y(ECon{x) 

Let us now finally introduce the main actor of this paper: 
Definition B.6 

set of the Pisot- Vijayaraghavan numbers: 

PV{A) := {.T e Az : X > 1 A (|y| < IVy G Con{x))} (B6) 
Let us recall that given a, 6 G K: 
Definition B.7 
floor of a: 

LOj := max{n G N : n < a} (B7) 

Definition B.8 

Modb(a) a-i_-jh (B8) 



Given a G M+ let us endow the interval [0, a) with the topology of the circle, id est by the topology induced by the 
following metric: 

Definition B.9 

-geodesic distance over [0, a); 

rf.eo.e.c(x,y):^(l^7l' , ;f|--yl^t; (B9) 

I a—\x — y\, otherwise. 

Given a sequence {a:„}„eN+ G [0,a)^+ and y G [0,a): 
Definition B.IO 

convergence modulo a: 

lim Xn = y {mod a) := Ve > 0,3N G N : idgeodesic{xn,y) < e Vn G N : ?i > TV) (BIO) 

n — *+oo 

Then [11], [zl: 
Proposition B.2 

Fundamental property of the Pisot- Vijayaraghavan numbers: 
HP: 



X G PV{A) (Bll) 

TH: 



lim Modi(x")=0 (B12) 

n— >+C50 
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PROOF: 

By the proposition lB.il 

+ if e Z yn e N+ : n > deg{x) (B13) 

y^Con(x) 

Modiix"- + y") = Vn e N+ : n > deg{x) (B14) 

y^Con{x) 

let us observe that since x G PV{A): 

V 2/" = O„^oo(-) (B15) 

y£Con{x) 

SO that: 

Modib" + Cl„^oo(-)] = Vn e N+ : n > deg(a;) (B16) 
n 

from which the thesis follows H 

Another important property of Pisot-Vijayaraghavan numbers concerns their link with suitable linear recurrence 
numeric sequences, according to the following: 

Definition B.ll 

linear recurrence numeric sequence: 

dif) 

{fnUeN e : (3 d{f) e N+, 3ci, • • • , c^if} e Z'^^-f^ ■ fn = Y ^''-^"-^ V?i G iV : ,i > dif)) (B17) 

k=l 

We will call the integer d(f) the degree of the linear recurrence numeric sequence f. 
Remark B.l 

It is important to remark that a linear recurrence numeric sequence is completely defined by primitive recursion once 
the values of /o, • • • , fd(f)-i have been assigned. 

Then (see the 13*'* chapter "Pisot sequences, Boyd sequences and linear recurrence" of (2^): 
Proposition B.3 

Linear recurrence numeric sequence associated to a Pisot-Vijayaraghavan number: 
HP: 



A e PV{K) (B18) 

TH: 



fn 

3{/ri}neN G linear recurrence numeric sequence : lim — — = A (B19) 



Remark B.2 
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Let us remark that, in in general, there might exist many primitive recursive numeric sequences associated to a 
Pisot-Vijayaraghavan number. 

Let us observe furthermore that: 
Proposition B.4 
HP: 



A e PV{A) : ord{X) = 2 
{/njnsN G 2^ hnear recurrence numeric sequence associated to A 

TH: 



3{c„}eM" : /„ = c„A" + 0„^+oo(-) A Um ^ = 1 (B20) 



Let y £ Con{\) : y 7^ A be the conjugate of A. 
Posed: 



PROOF: 



one obtains that: 

fn _ c„A" + 0„_>+oo(^) 



/„_i c„_iA"-i +0„^+oo(i) 
Using the fact that {/«} is associated to A: 



from which the thesis follows 



Remark B.3 



(B22) 



A + 0„^+oo(-) = ^A (B23) 

n Cn-l 



If the linear recurrence numeric sequence is of the form /„ = mfn-i + /n-2, with /o = and /i = 1 Binct has proved 
that the sequence c„ of the proposition IB. 41 is constant not only asymptotically: 

^ Vn e N (B24) 



Example B.l 

Given the golden number t := 

-Pt — minimal {-^^ — *^ 1 (B25) 

and hence Pr-minimai e 7'[Z]. 
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Furthermore r > 1 while: 

Con{T) = {--} (B26) 
r 

where | — i| < 1. 
Hence /e PV{A). 

The golden number is the least Pisot-Vijayaraghavan number of order two. 
Furthermore it is the least limit point of Pisot-Vijayaraghavan numbers: 

T = mmLP[PV{A)] (B27) 

A linear recurrence numeric sequence (of degree 2) associated to r is the sequence of Fibonacci numbers: 

Fn+2 ■= Fn+1 + Fn (B28) 

Fo := (B29) 
Fi := 1 (B30) 

Example B.2 

Given the plastic number p := ^ + (9+v^6 9)j^ . 

Fp—rriinirnali'^^ — ^ 1 (-^^-^) 

and hence Pp-mmimai € V[I^]. 
Furthermore p > 1 while: 

^ ,z(i + V3)(9- V69)W (-l-iV3)(9 + V69)3 {-I ^ i^/?>){% - ^/m)^ + i{i + ^/?>){% + Vm)^ . 

2-2^3§ 2-233§ 

(B32) 

^j^j^ , i{i+-/3){9-V69)i + (-l-i73) (9+769) i i _^ i (-l-»73)(9-769) S + ^(»+73)(9+769) ^ i ^ 

2-23 33 2-23 33 

Hence p e PF(A). 

The plastic number is the least Pisot-Vijayaraghavan number. 

Furthermore, introduced the notion of the height h{x) of an algebraic number x G A (a sort of measure of the 
algebraic complication needed to describe x for whose definition and properties we demand to the 3*'' chapter "Pisot 
and Salem Numbers" of [75|, to the section 3.5 "Smyth's theorem" of ^26(1 and to the 4*'* chapter "Small points" of 



[76|), the plastic number appears in the following: 

Theorem B.l 

Smith Theorem: 
HP: 



xGAz : x^O h- i Con{x) (B33) 

X 



TH: 
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A linear recurrence numeric sequence (of degree 3) associated to p is the sequence of the Padovan numbers: 

P(0) := P(l) := P(2) := 1 (B35) 

Pn Pn-2+Pn-3 (B36) 

Example B.3 

The number 1 + \/2 > 1 has as minimal polynomial: 

^i+v^-„.™/(^) = ^'~2x-l (B37) 

and is hence an algebraic integer. 
Furthermore: 



where obviously: 



Con{l + V2) = {I-V2} (B38) 



II-V2I < 1 (B39) 



Hence 1 + V2e PV{A). 

A linear recurrence numeric sequence (of degree 2) associated to 1 + V2 is the sequence of Pell numbers: 

ao := (B40) 

ai := 1 (B41) 

a„ := 2a„_i + a„_2 Vn G N : n > 2 (B42) 

We want to conclude this Number Theoretic section with a celebrated well known result, the recursive undecidability 
of Hilbert 10*'* problem showing its deep link with Theoretical Computer Science 0, 
Given a set C N: 

Definition B.12 

5* is Diophantine: 



3P e V[Z] : S = {xe'N : P{x) = 0} (B43) 



Tiieorem B.2 

Matiyasevich Theorem: 
HP: 



S" C N (B44) 



TH: 



S is Diophantine if and only if it is recursively enumerable 

Theorem IB.2I implies that it doesn't exist an algorithm that, receiving as input a polynomial P £ V{'&) outputs 
"yes" if its has non-negative integer solutions while it outputs "no" if it doesn't. 
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APPENDIX C: A BRIEF INFORMATION THEORETIC ANALYSIS OF SINGULAR 
LEBESGUE-STIELTJES MEASURES SUPPORTED ON CANTOR SETS AND ALMOST PERIODICITY 

Demanding to [t^I, [ll], [l^ for any further information let us recall that given a continuous function / : M i-^ R 
non decreasing, id est such that: 

I{x)>f{v) Vx,2/eM : x>y (CI) 

we can introduce the following: 
Definition C.l 

Lebsegue-Stieltjes measure associated to the f: 
the measure /// : B{M.) i— > [0, +00) such that: 

fif[{a,b)] lim /(6-e) - lim /(a + e) Va, 6 G K : a < 5 (C2) 

where B{M.) is the Borel-cr-algebra on M, id est the cr-algebra induced by the topology induced by the usual metric 
d{x,y) := |a; - y\. 

For instance, choosing the identity function /d : M ^ M: 

Id{x) := X (C3) 

the Lebesgue-Stieltjes measure fiid is nothing but the Lebesgue measure on numbers VA{^J.Lebesgue) (where ^Lebesgue 
is the Lebesgue measure on sequences of the definition lA.Sp . 
Given a measure /i : S — > [0, +00) let us recall that: 

Definition C.2 

set of pure points of ^ : 

PP{y) := {a; e M : ^J.{{x}) ^ 0} (C4) 

Definition C.3 

/i is pure point: 

fi{X) = ^ /i({x}) \fX e B{R) (C5) 

xex 

Definition C.4 

fi is absolutely continuous (with respect to the Lebesgue measure): 

H -< fiM ■■= KB) = ^ hm{B) = (C6) 

Let us recall the basic: 

Tileorem C.l 

Radon-Nikodym Theorem: 
HP: 



TH: 



There exists a /i/^-measurable function, called the Radon-Nikodym derivative of with respect to ///d and 
^,suchthat4/(a;)dM^) = /k/(x)^ 



consequentially denoted as -j^, such that Jjg f{x)dfi{x) = f{x)-^p^{x)dfiid{x) for every /x-measurable function f 
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Definition C.5 

fi is singular (with respect to the Lehesgue measure): 

M ± Hid 35 : ti{S) = A - S) ^ Q (C7) 

Let us then recall the following: 
Theorem C.2 

Lehesgue Decomposition Theorem: 
HP: 



H : B{M) ^ [0, +oo) measure 

TH: 



3\{Hp.p, Ha.c, fJ'sing) triple of measures : /ip.p is pure point A ^ ^ ^a.c. A ^ ± fising 
Given a map / : [a, b] M: 
Definition C.6 

/ is absolutely continuous: 

n 

Ve>0, 35>0:^|/(6,)-/K)| <e 

n 

V{(a»,fe»)}r=i : (az,60 C (a,6) A (a„6,)n (aj,6j) =0 ^IId[{a^,b,)] < 5 \fn e N+ (C8) 

1=1 

A remarkable feature of Lebesgue-Stieltjes measures is the following: 
Theorem C.3 

On the singularity of Lebesgue-Stieltjes measures: 
HP: 

Hf Lebesgue-Stieltjes measure 

TH: 

/i/ is an absolutely continuous measure if and only if f is an absolutely continuous function 
Let us now show how Algorithmic Information Theory enters in the game. 

Given a trinary alphabet A := {oi, 02, 03} let us consider the set of all sequences over A not containing the letter 
02: 

NOA,a, {x e A'*+ : \x\a, = 0} (C9) 
As we have already observed in the remark |IH.4[ the theorem IHL4I implies that: 

NOam2^RANDOM{A^+) = (CIO) 
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Let us now consider the binary alphabet B := {ai, a^}. 
Clearly: 

NOA,a2 = ^ RANDOM{B^^+) (CU) 

At a first sight the eg uat ion I C 1 01 and the equation [CTl] might appear incompatible. 

A deeper analysis allows, anyway, to understand that everything is absolutely consistent: 

simply we have to remember that the correct way of passing from the alphabet A to the alphabet B is through the 
transition map Ta.s of the definition lA. 61 
Actually the theorem IIII.3I guarantees that: 

RANDOM{B^^+) = Ta.b[RANDOM{A^^+)] (C12) 



Remark C.l 

In a very interesting paper [sO] that has greatly inspired this section, Cristian Calude, Ludwig Staiger and Karl Svozil 
have introduced a generalization of the notion of algorithmic randomness, id est of the definition IIII.241 to sequences 
{xn} each letter of which belongs to a different alphabet A„. The analysis of the probably deep link existing between 
such a generalization and what we are going to discuss here is something that, at present, we are not able to formalize 
and is left as an interesting challenge for the reader. 

Let us now introduce the set: 
Definition C.7 

Cantor's middle third set: 



CA^a, ■■= VAiNOA.a,) (C13) 



It may be easily proved that: 
Proposition C.l 
1. 

2. 

3. 



t^IdiCA^a,) - (C14) 



\CA.aJ > Ho (C15) 



dimHausdorff{CA,a2) = , ^ I ■ ! (C16) 

log \A\ 

Given n E N+ let us introduce the following: 

Definition C.8 

Cantor function: 
the map /a, a : B^ i-^ R: 

fAA^) ■■= vb{x) + J^Vf e B+ : vb{x) E {va(x) + jj^,va{x) + (C17) 

It may be easily proved that the map / : [0, 1] i-^ R: 

/>) firBix)) (C18) 



is a continuous non-decreasing function so that it induces a Lebesgue-Stieltjes measure /i^. 
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Since / is not absolutely continuous it follows that fij is singular with respect to the Lebesgue measure; 
actually: 

fif[CA,a,] = 1 (C19) 

while we know by the proposition IC . 1 1 that CA,a2 has vanishing Lebesgue measure. 
Let us remark that C^.aa is a Cantor set in the sense of the definition lA. 21 

Now the measure /xj : B[[0, 1)] i-^ [0, +oo) induces naturally the measure fif : ;B[y4^+] i— > [0, +oo) defined by: 

HfiX) := fijivAiX)) (C20) 

Then: 

fJ-Lebesgue,B[NOA,a2] = ^^Lebesgue,B[B^*] ~ 1 (C21) 
[iLehesgu.A^OA.a2^RANDOM{B^'^)\ = ^iebe.gne.B [i?A^I?OM )] = 1 (C22) 

[A^OA,aJ = (C23) 
p.j\NOA,a.^ = 1 (C24) 

where we have used the theorem lA.il 

The whole story may be easily generalized in the following way: 

given a finite alphabet A : 2 < \A\ < Hq and a letter a e yl let us introduce the set of all the sequences on A not 
containing the digit a: 

NOA,a {x e : \x\a = 0} (C25) 
and the new alphabet B obtained from A excluding the letter a: 

B ■= A-{a\ (C26) 

Then one can introduce the set: 
Definition C.9 

generalized Cantor set associated to {A, a): 

CA,a ■■= VA{NOA,a) (C27) 

Given n e N+ let us introduce the following: 
Definition C.IO 

generalized Cantor function associated to {A, a): 
the map fA,a '■ ^ ^■ 



\B\ „^ , 1 \B\ 



(C28) 



IaA^) ■■= VBix) + j^^yx e B+ : VBix)eivAix) + jj^,VAix) + jj^) 
and the associated map / : [0, 1] i-^ M: 

IaAx) := f{rB{x)) (C29) 

Then: 
Proposition C.2 

Properties of the generalized Cantor functions: 
HP: 
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< lex{a) <\A\-l 



(C30) 



TH: 



1- /a, a is a continuous non decreasing function 

2- /a, a is not absolutely continuous 

PROOF: 

It is sufficient to follow step by step the analogous proof holding in the particular case of the Cantor middle-third set 
and generalize in the natural way H 

The proposition IC . 21 implies that one can introduce the Lebesgue-Stieltjes measure /ij. Then: 

Corollary C.l 

is singular 
PROOF: 

The thesis follows by the proposition IC.3I and the proposition IC.2I ■ 
Furthermore: 



fJ-Lebesgue,B[NOA,a] = fJ-Lebesgue,B[B +] — 1 (C31) 

fiLebesgue,BWOA,anRANDOMiB^^)] = f^LebesgueMR^N DOM (3^+)] = 1 (C32) 

tJ-Lebesgue,A[NOA,a] = (C33) 

M[NOA,a] = 1 (C34) 



We will now briefly outline the deep link existing between the singular Lebesgue-Stieltjes measures associated to 
generalized Cantor functions and the theory of almost periodic functions, sequences and measures [8i| . [s^ . [83| . 
Given a map / : R C: 

Definition C.ll 

/ is periodic: 

3rG(0,+oo) : f(t + T) = f(t) yteR (C35) 

Given a periodic function f: 
Definition C.l 2 

fundamental period of f: 

Periodfund[f] ■■= min{r e (0, +00) : f{t + T) = f{t) Vi G M} (C36) 
We will denote by PERt{R, C) the set of all the periodic functions. 
Remark C.2 
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In the section Hill we saw that the combinatorial information and the algorithmic information of a periodic sequence 
(according to the definition IIII.7P are very low since a periodic sequence is completely specified assigning the values 
it takes over a single period. 

The same thing can clearly be said of periodic functions. 

Clearly: 
Proposition C.3 

PERt{M.., C) is a complex Hnear space VT £ (0, +00) (C37) 

Given n G N+, uji, ■ ■ ■ ,LUn G (0, +00) and a map / : M i-^ C: 
Definition C.13 

/ is quasiperiodic with pulsations lui, ■ ■ ■ ,LUn: 

{34> : E" 1-^ C : (j){xi, ■ ■■ ,Xi + 2TT,-- ■ ,Xn) = 4>ixi, • • • , a;^, • • • , a;„) V« G {1, • • • , rt}) A 

(fit) = (/.(a;it, 

Definition C.14 

uJi, ■ ■ ■ ,ujn are rationally independent: 

n 

{J2 hu^i = => ki = --- = kn = 0) V/ci, • • • , fc„ e Q 
1=1 

Then: 
Proposition C.4 

On the relation between periodicity and quasiperiodicity: 

1. f is periodic f is quasiperiodic 

2. f is periodic <^ f is quasiperiodic with rationally dependent pulsations 
Definition C.15 

f is a complex trigonometric polynomial: 

n 

3n e N+, , • • • , A„ e C, 3wi, • • • , w„ e (0, +c») : f{t) = ^ ^fc exp{iuJkt) Vt G M (C40) 

k=l 

Proposition C.5 

/ complex trigonometric polynomial / is quasi-periodic (C41) 

PROOF: 

Given the complex trigonometric polynomial: 

n 

fit) := ^^feexp(iwfet) (C42) 
fe=i 

let us introduce the following map ^ : E" 1— > C: 

n 

(j){xi,- ■ ■ ,Xn) := ^Akexp{ixk) (C43) 

k=l 



■■■ ,ujnt) VteM) (C38) 



(C39) 
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Obviously: 

(j>{xi + 27r, • • • ,x„) = • • • = • • • ,a;„ + 27r) Vxi, • • • ,x„ G M (C44) 

and: 

f(t) = (/.(wii,- • • ,u;„i) VieM (C45) 

■ 

Let us now introduce Harald Bohr's notion of almost periodic function. 

Definition C.16 

almost periodic function: 
a map / : M i-^ C : 

Ve > 03Tc complex trigonometric polynomial : \f{x) - T^{x)\ < e V.t G R (C46) 
We will denote hy A — PER{R, C) the set of all the almost periodic functions. 

To appreciate the meaning of the definition IC . 161 let us recall that [77] given a metric space (M, d) and a set S C M: 

Definition C.17 

S is dense in {M,d): 

yyeS 3{xn}neN G M"* : Hm d(a;„, y) = (C47) 

n — >+oo 

Then definition I C . 1 61 implies that: 
Proposition C.6 

{ complex trigonometric polynomials } is dense in the metric space (A — PER{R, C), di) 

where: 

:= sup 1/(0 -ffWI (C48) 

teR 

Furthermore: 
Proposition C.7 

On the relation between quasiperiodicy and almost periodicity: 

1. f is quasiperiodic =^ f is almost periodic 

2. f is almost periodic ^ f is quasiperiodic 

The basic properties of almost periodic functions are encoded in the following: 
Proposition C.8 

Basic properties of almost periodic functions: 
1. 

/ is almost periodic <^ Ve G (0, +00) 3/(e) e (0, +00) : 

(VaeM,3eeR : \f{x + £.)-f{x)\ < eVx € (a, a + /(e)) (C49) 



2. f complex trigonometric polynomial =4> f almost periodic 

3. A — PER{R, C) is a complex linear space 



Km f{x)dx = (C51) 

da T^+oo . 
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4. 

/i • /2 is almost periodic V/i, /2 almost periodic (C50) 

5. 

d 

Given an almost periodic function f: 
Definition C.18 

mean value of f: 

i-a+T 

M[J\ lim / f{x)dx (C52) 

T^+ca J ^ 

Then: 
Proposition C.9 

Basic property of the mean value of an almost periodic function: 

supportujM[e-x:p{~'iLLix)f{x)\ < Nq (C53) 

Given a generic function f <C: 
Definition C.19 

Fourier transform of f: 

_ j limT^+oc Jl^^ dweyi-p{-iujx)f{x), if the limit exists; fnKA\ 
I +O0, otherwise 

The proposition IC.9I implies that: 
Proposition C.IO 

Fourier spectrum of almost periodic functions: 
HP: 



f almost periodic function 

TH: 



1. 

+CSO 

3{A„}eC^3{cc>4e(0,+oof : = ^ A„5(co-co„) (C55) 

n— — oo 

2. if f is periodic of fundamental period ^ then lj„ = nui Vn G Z 

PROOF: 

1. The thesis immediately follows combining the definition IC.19I and the proposition IC. 91 

2. The thesis is nothing but the basic theorem concerning Fourier series of periodic functions 
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As to the existence of the Fourier transform of almost periodic functions there exists a very powerful result [83 | : 
Proposition C.ll 

About the specification of an almost periodic function through its Fourier transform: 
HP: 



M e (o,+oo)^ 

+00 
n— — oo 

TH: 

-\-OQ 



3 f e A - PER{R,C) : = ^ A„(5(tj - w„) 



Let us now consider a generic measure /i : S(M) [0, +00) on the real line and let us introduce the following: 
Definition C.20 



Fourier transform of fi: 



Obviously: 
Proposition C.12 
HP: 



TH: 



^[/i](u;) := / e-*""d^(x) (C57) 



•^M = n:r^] (C58) 

where the Fourier transform of the left hand side is intended in the sense of the definition IC.20I while the Fourier 
transform of the right hand side is intended in the sense of the definition IC.19I 

PROOF: 

It is sufficient to observe that: 

= / e~^'---^{x)d^lM{x) - e-''^'- 4^{x)dx = ^[-^] (C59) 

jR a^iid J_oo a^iid anid 

m 

The situation is, instead, more complicated in the case in which the measure fi is singular with respect to the 
Lebesgue measure d/iid- 

Let us consider in particular the case in which — l^fj^ is the singular measure associated to a generalized Cantor 
function /a,q. 

Though we have not yet a proof, we think that Eberlein's Theorem (see the 11*'* chapter "Almost periodicity of 
the generalized Fourier transform " of [s^) strongly supports the following: 

Conjecture C.l 

About the Fourier transform of the Lebesgue- Stieltjes measure associated to a generalized Cantor function: 

■^[^/a a\ almost periodic (C60) 
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APPENDIX D: QUASICRYSTALS 

Given n e N+ let us introduce the following: 

Definition D.l 

n- dimensional euclidean space: 

the Hilbert space E" := (M", •), where x ■ y :— J27=i ^iVi is the usual euchdean inner product. 
Given E" and r e M+: 

Definition D.2 

ball with center x and radius r: 

Br{x) := {yeE" : |y-f| <r} (Dl) 
Given S* C E" let us introduce the following f^: 
Definition D.3 
S is a Delone set in E" : 
1. 



2. 



3ro e M+ : \x - y\ > 2ro yx,y e S x ^ y 



3Ro e M+ : Br{x) n 5 7^ Vr > i?o , Vf G E" (D2) 



Following the generalization of the mathematical definition of a crystal given in [ISli [2^ the score of 
the indications given in 1992 by the Commission on Nonperiodic Crystals established by the International Unit of 
Crystallography, we will characterize a crystal by the condition that its diffraction pattern exhibits a countable infinity 
of Bragg peaks. 

Given a Delone set S: 

Definition D.4 

density distribution of S: 

Ps{x) := ^5{x-s) 

Definition D.5 

autocorrelation of S: 



jsix):^ Pix)*pi-x) (D3) 

where ★ denotes the convolution operator: 

{f*g){x) := f dyf{x)g{x-y) (D4) 



Definition D.6 

spectrum of S: 

the measure /is on (E",S(]R")) identified by limL^+00 (2L)" -^fen[-L,L]"(3^)] 
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where B{W^) is the Borel-cr-algebra of R" and where T denotes the Fourier transform [73|, |84| . 

By the Lebesgue Decomposition Theorem (i.e. the straigthforward multidimensional generalization of the theorem 
IC.2[) we know that /xg may be decomposed uniquely as: 

where fJ-^g'^'"^ is a pure point measure, /i^"'^''' is absolutely continuous w.r.t. the Lebesgue measure fJ-Lebesgue while 
^Oi^s) jg singular w.r.t. fJ-Lebesgue- 

The mathematical characterization of crystals involves only /i^'^'^; introduced the following [tt} : 
Definition D.7 
sef o/ pure points of /ig : 

PP(/is) {^eK" : /.s({f})^0} = {f G R" : M^' ^-^ax}) «} (D6) 

we can at last define: 
Deflnition D.8 

S is a crystal: 

\PPi^^s)\ = Ho (D7) 

Remark D.l 

It is claimed in the book review [H of 01 that the definition |Dj is not mathematically rigorous. 

It should be noticed, with this regard, that the author of such a criticism hasn't proposed [86j an, according to him, 
better definition. 

Demanding to 2J| for further details let us observe that since for every Delone set S: 

-^nisn[-LM^)] e 5'(M") Vi e M+ (D8) 

its limit for L — > oo is perfectly well defined with respect to the natural topology of the Schwartz space of tempered 
distributions (tt} . 

As to definition ID. 81 it is important to remind the double nature of the Dirac delta as a tempered distribution and 
as a Lebesgue-Stieltjes measure. 

Let us recall that given k £ N+ vectors ai , • • • , e M" : 

Definition D.9 

li-module generated by ai, ■ ■ ■ , 0^; 

fc 

r(ai,---,afe) := {^niidi nii, e Z i ^ 1, ■ ■ ■ , k} (D9) 

i=l 

Definition D.IO 

lattice in E" .■ 

a Z-module generated by a basis of E" 

Let us now consider the particular case in which S is a lattice T{di, • • • , a„) (where hence {oi, • • • , a„} is a basis of 
E") and let us introduce the following: 

Definition D.ll 

dual lattice ofT{di, • • ■ ,dn)'- 

r'(ai, • • • ,a„) := T[d\, ■■■ ,a^) 
where {a!^, ■ ■ ■ , aj^} is called the basis dual to the basis {ai, • • • , a„} and is defined by the condition: 



Oj ■ a 



■' := (5,,, = - ,n (DIO) 



A corner stone of Mathematical Crystallography (both classical and quasi) is the following Q : 
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Theorem D.l 

Poisson's summation formula: 

•?^(Pr(ai,---,a„)) = Pr'(ai,-,a„) (Dll) 

from which it foUows that: 
Corollary D.l 

r(ai, • • • , a„) is a crystal 

The problem of characterizing which nonperiodic Dclone sets are crystals is an extremely high and steep mountain 
that has been challenged by some of the best minds in the scientific community (see e.g. 0] and [23l|). 

Let us now introduce the following: 

Definition D.12 

isometry in E" ; 

a map : E" E" such that: 

\cj){S-y)\ - \x-y\ Vx,yeE" 

The isometrics of E" form a group that we will denote by /s(E"). 
Given a Delone set S: 

Definition D.13 

symmetry group of S: 

Sim{S) := {0 G /s(E") : (j) o = t^s} (D12) 
Given a group G and an element g ^ G: 
Definition D.14 
order of g: 

ord{g) := min{?i e N+ : 5" = 1} 



Definition D.15 

Hiller's function: 

Hil-.N+i-^N : Hilik) := min{n e N : .g G GL{n, Z) A ord{g) = k) 
Introduced the well-known: 
Definition D.16 

Euler's <j) function: 

: N+ N+ : (t>{n) := \{k en+ : k < n A gcd{k, n) = 1}\ 

then [131 : 
Theorem D.2 
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Hiller's Theorem: 

oo 

n=l {n6N+ :p°"5<t2} 

where p„ is the n^^ prime number. 

The Hiller function is algorithmically implemented in the section [E) 
It may be easily verified that: 



n 


Hil{n) 


1 





2 





3 


2 


4 


2 


5 


4 


6 


2 


7 


6 


8 


4 


9 


6 


10 


4 


11 


10 


12 


4 


13 


12 


14 


6 


15 


6 


16 


8 


17 


16 


18 


6 


19 


18 


20 


6 


21 


8 


22 


10 


23 


22 


24 


6 


25 


20 


26 


12 


27 


18 


28 


8 


29 


28 


30 


6 


31 


30 


32 


16 


33 


12 


34 


16 


35 


10 


36 


8 



Remark D.2 

Let us remark that the constraint p^" 7^ 2 in theorem ID. 21 is essential as it is proved in where, contrary to 0], 
such a theorem is reported correctly. 

It can be easily verified that: 



67 



Corollary D.2 

Hil{n) < 3 ^ n e {1,2,3,4,6} 

Given a Delone set S in E" : 
Proposition D.l 

Crystallographic restriction: 

S* is a lattice ^ Hil{ord{g)) < n Vg G Sim{S) 
PROOF: 

The thesis follows by the same definition ID. 151 of Hiller's function. ■ 
Definition D.17 

S is a quasicrystal: 

S* is a crystal A 3^ e Sim{S) : Hil{ord{g)) > n 

Example D.l 

By the coroUarv ID . 2l n-fold rotational symmetry axes, i.e. axes around which there is rotational symmetry for rotations 
by can exist in two and three dimensions if and only if n G {2, 3, 4, 6}. Hence a bidimensional or three-dimensional 
crystal exhibiting an n-fold rotational symmetry axis with n ^ {2, 3, 4, 6} is a quasicrystal. 

Example D.2 

Passing from three-dimension to four-dimension the only symmetries that become allowed are those of order 5,8,10,12. 
Hence a four-dimensional crystal exhibiting a symmetry whose order ^ {2, 3, 4, 5, 6, 8, 10, 12} is a quasicrystal. 
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APPENDIX E: MATHEMATICA IMPLEMENTATION OF THIS PAPER 

The notions introduced in this paper may be implemented algorithmicaUy through the following Mathematica 5 
[H, [iil notebook: 

Off [General : : spelll] ; 

$MaxExtraPrecision=\ [Infinity] ; 

<<DiscreteMath'Combinatorica' ; 

<<NumberTheory ' AlgebraicNumberFields ' 

(***words [alphabet , 

n] is the list of all the words of length n in lexicographic ordering***) 

words [alphabet_ ,n_] : =Strings [alphabet ,n] 

(***wordsupto [alphabet , 

n] is the list of all the words of length less or equal to n in \ 
lexicographic ordering***) 

wordsupto [alphabet_,n_] : = 

If [n\ [Equal] 1 , words [alphabet , 1] , 

Join [wordsupto [alphabet ,n-l] , words [alphabet ,n] ] ] 

(*** subwords [x] gives the list of all the sub- 
words of the word x in lexicographic ordering***) 

subwords [x_] : =Flatten [Table [Tcike [x , {i , j }] , {i , 1 , Length [x] } , { j , i , Length [x] }] , 1] 

(^^^ generalizedselect [1, predicate, 

y] picks out all elements e of the list 1 for which the binary \ 
predicate predicate [e ,y] is True ***) 

(*** WARNING: 

since the instruction is implemented throwing away elements of pattern \ 
String the list 1 must not to contain strings ***) 

generalizedselect [1_ ,predicate_ ,y_] : = 
DeleteCases [ 

Table [If [predicate [Part [l,i] ,y] ,Part[l,i] ."throw away"] , {i , 1 .Length [1] >] . 
x_String] 

havinggivenlengthQ [word_.n_] :=Equal [Length [word] .n] 

(*** language [word, n] is the language of length n of the word x 
***) 

Isinguage [word_.n_] : =generalizedselect [subwords [word] .havinggivenlengthQ .n] 

(*** topologicalentropy [alphabet .word. 

n] is the approximation at level n of the topological entropy of a \ 
sequence ***) 

topologicalentropy [alphabet_ . word_ . n_] : = 

Divide [Log [Length [alphabet] . Length [language [word. n] ] ] .n] 
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occurences [x_,y_] :=If [AtomQ[x] , Count [{x} , y] .Count [x,y]] 

(*** substitution [rule, 

x] is the action of the substitution rule over the word x ***) 

substitution [rule_ , x_] : =Flatten [Map [rule , x] ] 

incidencematrix [alphabet. , rule_] : = 

Table [occurences [rule [Part [alphabet , j] ] ,Part [alphabet , i] ] ,{i, 1 , 
Length [alphabet] }■ , { j , 1 , Length [alphabet] }•] 

(*** sequence [rule, x,n] is the n- 

the iterate of the action of the substitution rule over the word x ***) 

sequence [rule_ , x_ , n_] : = 

If [n\ [Equal] 1 , substitution [rule, x] , substitution [rule , sequence [rule ,x,n-l] ] ] 

f ibonacci [x_] : =If [x\ [Equal] , {0 , 1} , 0] 

padovan [x_] : =If [x\ [Equal] , {1 , 2} , If [x\ [Equal] 1,2,0]] 

padovannumber [n_] : =If [n\ [LessEqual] 2,1, padovannumber [n-2] +padovannumber [n-3] ] 

procedure [rule_,x_,n_] :=Do [Print [sequence [rule ,x,k] ] ,{k,l,n}] 

algebraicorder [x_] := 

If [x\ [Element] Algebraics .Exponent [MinimalPolynomial [x,y] ,y] ,\ [Infinity] ] 

dif f erentQ [x_,y_] : =Unequal [FullSimplif y [x-y] .0] 

conjugates [x_] : = 

If [x\ [Element] Algebraics . 
generalizedselect [ 

Table [Root [MinimalPolynomial [x] ,k] ,{k. 1 .algebraicorder [x] }] .dif f erentQ. 
x] ,\ [Infinity]] 

absolutevaluelessthanoneQ [x_] :=Abs [x] <1 

pvnumberQ [x_] : = 

And[AlgebraicIntegerQ [x] ,x>l , 

Length [Select [conjugates [x] , absolutevaluelessthanoneQ] ] \ [Equal] 
Length [conjugates [x] ] ] 

polynomial [x_,listof coff icients_] := 

Sum [listof coefficients [ [i]] *x"i,{i, 1, Length [listof cof f icients] }] 

\ [Tau] =GoldenRatio ; 

\[Rho]=Root[(#-3-#-l)&,l] ; 

pointofpolygon[k_, i_] := {Cos [(2*\ [Pi] *i) \k] , Sin[(2*\ [Pi] *i)/k] } 



lineof polygon [k_, i_, j_] : =Line [{pointof polygon [k, i] ,pointof polygon [k,j] }] 

polygonpent agram [k_] : =Flatten [Table [lineof polygon [k , i , j ] 
,{i,l,k>,{j,l,k}]] 
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polygonpicture [k_] : =Show [Graphics [polygonpentagram [k] ] ] 

poiiitofpentagrEmi[n_, i_] : = 

If [n \ [Equal] 1, {Cos [(2*\ [Pi] *i) /5] , 

Sin[(2*\[Pi]*i)/5]>, - (GoldenRatio/(l + 2*GoldenRatio) ) 
pointof pentagram [n - 1, i]] 

lineof pent agram [n_ , i_ , j _] : =Line [{pointof pentagram [n , i] , pointof pentagram [n , j ] }] 

pentagram [n_] : =Flatten [Table [lineof pentagram [n, i , j ] 
,{i,l,5},{j,l,5}]] 

pentagramipicture [n_] :=Show [Graphics [Table [pentagram [i] ,{i , 1 ,n)-] ] ] 

pvpoint [\ [Lcunbda] _ ,n_] : = 
If [n\ [Equal] 1,{1, 

0> , {Re [Exp [\ [Imaginaryl] *Mod [ (2*\[Pi]) * (\ [Lambda] "n) , 2\ [Pi] ]]] , 
Im [Exp [\ [Imaginaryl] *Mod [ (2*\ [Pi] ) * (\ [Lambda] "n) , 2\ [Pi] ] ] ] }] 

pvside [\ [Lambda]_,n_] :=Line [{pvpoint [\ [Lambda] ,n] , pvpoint [\ [Lambda] ,n+l]}] 

pvline [\ [Lambda] _,i_,j_] :=Line [{pvpoint [\ [Lambda] ,i] .pvpoint [\ [Lambda] ,j]}] 

pvcurve [\ [Lambda] _ ,n_] :=Table [pvside [\ [Lambda] ,k] ,{k,l,n-l}] 

pvcurvepicture [\ [Lambda] _ , n_] : = 

Show [Graphics [{Circle [{0,0} , 1] , pvcurve [\ [Lambda] ,n]}]] 

pvpentagreim [\ [Lambda] _ ,n_] : = 

Flatten [Table [pvline [\ [Lambda] ,i, j] ,{i,l,n},{j ,l,n}]] 

pvpentagreimpicture [\ [Lambda] _ ,n_] : = 

Show [Graphics [{Circle [{0,0>, 1] , pvpent agr am [\ [Lambda] ,n]}]] 

elementaryrotation [alphabet_ , spacings_ , letter_ , initialpoint_] : = 

Mod [Part [spacings , Extract [Position [alphabet , letter] ,{1,1} ] ] +initialpoint , 
2\ [Pi]] 

orbit [alphabet_ , spacings_ ,word_ ,initialpoint_,n_] := 
If [n>Length [word] , "undefined" , 
If [n\ [Equal] 1 , 

elementaryrotation [alphabet , spacings ,Part [word, 1] , initialpoint] , 
elementaryrotation [alphabet, spacings, Part [word,n] , 
orbit [alphabet , spacings , word , initialpoint ,n-l] ] ] ] 

pointof pisotsequence [alphabet_ , spacings_ ,rule_,x_ , 
n_] : ={Re [Exp [\ [Imaginaryl] * 

orbit [alphabet , spacings , sequence [rule,x,n] ,0, 
Length [sequence [rule , x , n] ] ] ] ] , 
Im [Exp [\ [Imaginaryl] * 

orbit [alphabet, spacings, sequence[rule,x,n] ,0, 
Length [sequence [rule ,x,n] ] ] ] ] } 

curveoff pisotsequence [alphabet. , spacings_ , rule_ , x_ , n_] : = 
Flatten [Table [ 

Line [{pointof pisotsequence [alphabet , spacings , rule ,x,k] , 

pointof pisotsequence [alphabet , spacings , rule ,x,k+l] }] , {k, 1 ,n-l}] ] 
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pisotsequencepicture [alphabet _ , spacings_,rule_,x_,n_] := 
Show [Graphics [{Circle [{0,0>,1] , 

curveoff pi sot sequence [alphabet , spacings , rule , x ,n] }] ] 

Hiller[n_] := 
If [n==l,0, 
If [n==2,0, 

If [Length [Factor Integer [n] ] ==1 , 

EulerPhi [Factor Integer [n] [ [1] ] [ [1] ] "Factorlnteger [n] [ [1] ] [ [2] ] ] , 
Sum[Hiller[ 

Factorlnteger [n] [ [i] ] [ [1] ] "Factorlnteger [n] [ [i] ] [ [2] ] ] , {i , 1 , 
Length [Factorlnteger [n] ] }] ] ] ] 
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